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ABSTRACT 

We assess the metal content of the cool (^10^ K) chxumgalactic medium (COM) about galaxies at z < 
1 using an H I-selected sample of 28 Lyman limit systems (LLS, defined here as absorbers with 16.2 < 
logTVni ^ 18.5) observed in absorption against background QSOs by the Cosmic Origins Spectrograph 
on-board the Hubble Space Telescope. The A'hi selection avoids metallicity biases inherent in many 
previous studies of the low-redshift CGM. We compare the column densities of weakly ionized metal 
species (e.g., O II, Si II, Mg II) to A^hi in the strongest H I component of each absorber. We find that 
the metallicity distribution of the LLS (and hence the cool CGM) is bimodal with metal-poor and 
metal-rich branches peaking at [X/H] ~ —1.6 and —0.3 (or about 2.5% and 50% solar metallicities). 
The cool CGM probed by these LLS is predominantly ionized. The metal-rich branch of the population 
likely traces winds, recycled outflows, and tidally stripped gas; the metal-poor branch has properties 
consistent with cold accretion streams thought to be a major source of fresh gas for star forming 
galaxies. Both branches have a nearly equal number of absorbers. Our results thus demonstrate there 
is a significant mass of previously-undiscovered cold metal-poor gas and confirm the presence of metal 
enriched gas in the CGM of z < 1 galaxies. 

Keywords: cosmology: observations — galaxies: halos — galaxies: abundances — - galaxies: kinemat- 
ics and dynamics 



1. INTRODUCTION 

One of the most pressing problems in galaxy forma- 
tion is to understand how gas accretion and feedback 
influence the evolution of galaxies and the intergalactic 
medium (IGM). Modern cosmological simulations can- 
not explain the obser ved mass-metallicity relationship 
(Tremonti et al. 2004) and color bimodalities observed 
in galaxies withou t invoking such large-scale flows (e.g., 
Keres et al. 2005). Stars cannot continue to form in 
galaxies over billions of years without a replenishment of 
gas from the IG M (e.g., |Maller fc Bullock[2004[ |Dekel fc 
Birnboim|2006 ), while feedback from star tormat ion and 
AGIN activity can fuel massive outflows from a gas-rich 
galaxy that may choke off star formation (e.g. , |Kacprzak 
et al.||20081 J Oppenhei mer et al.|[20T0l [20T2| l EumagaUi 
et al.pOlTb Faucher-Giguere et al.||2011[ ). The compe- 
tition between these regulating processes is played out 
in the circumgalactic medium (CGM), the interface be- 
tween galaxies and the IGM, where outflows may cool 
and stall, and where infalling gas may encounter winds 
or be shock-heated and halted. 
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Cosmological simulations have indicated the impor- 
tance of feedback a nd accretion through the CGM for 
several decades (e.g.. White fc Frenk|1991 Dekel & Birn- 



boim 2006). Modeling these phenoniena with more real- 
istic physics and modified methods to solve the hydrody- 



computational power (e.g.. 


Springel & Hcrnquist 2003 
ler et al..2010, Springel.20l'n 


Kobayashi|2004 lOppenheinr 


Smith et al.||2011 


Keres et al.||2012 |Sijacki et al.|2012 


Vogelsberger et a 


2012). All these simulations show in 



Morris et al. 



1991 
I'ripp et^ 



excellent tracers of the CGM of galaxies within ^300 kpc 
(about the virial radius, jFord et al.^2013) . Several obser- 
vational results from the analyses of the galaxy- absorber 
connection also show that the virialized CGM is traced 
by the strong Lyg absorbers (e.g., 
geron et al.||199 4: Lanzctta ct al 
1998; Pcnton et al.ll'2002; Bowen ct al. 2002; Chen et al 
2005; Morris & Jannuzi 2006; Wak ker fc Sav age 2009r 
jProchaska et al.„2011b) . These observations show that 
there is significant clustering between galaxies and strong 
H I (TVhi ^ 10^^ cm~^) absorbers, which points to a 



In contrast, the weak H I 
'^) are at best very loosely 



strong physical connection, 
absorbers (A^hi ^ 10^'* cm 
connected to galaxies, as evidenced by a weak or absent 
clustering signal in the galaxy-absorber two-point cross- 
correlation function. 

The most sensitive way to study the large-scale CGM 
around galaxies is therefore to investigate the properties 
of the strong H I absorption (and associated metal lines 
associate d with it) along QSO sigh tlines, as first pro- 
pos ed by Bahcall fc Spitzer ( 1969[ ) and then observed 
by Bergeron ~( 1986 r Here, we focus on the Lyman 



limit systems (LLS) as tracers of the C GM (e.g., Tytler 
1982[ |Steidel|1990"l [Bergeron et al.|1994[ ) , absorbers with 
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Lehner et al. 



large enough surface density of H I atoms to produce a 
decrement in the UV flux at the Lyman hmit. Tradi- 
tionaUy, they have been separately categorized as par- 
tial Lyman systems (16.2 < logTVni < 17) and LLS 
(17 < logiVni < 19). In this paper, we will just de- 
fine both categories as LLS for simplicity. Since all these 
absorbers trace CGM gas, in this context there is not a 
need to distinguish absorbers with tll < 1 (partial-LLS) 
and Tll > 1 (LLS), where tll is the optical depth at 
the Lyman limit. The importance of the LLS for char- 
acterizing the CGM and phenomena therein (large-scale 
flows) has been demonstrated in recent cosniol ogical sim- 
ulations (e.g.,[Faucher-Giguere fc Kere§r20lTbl|FumagaUi 
[St 



et al.|20lTbl [Stewart et ai.|2Ulla[|van de Voort fc Schaye 



2012r|Shen et al. p013p as well as observationally as we 
now discuss in more detail. 

Two recent results from our group have not only con- 
firmed that LLS probe the CGM, but also demonstrated 
they can trace wi dely varying proce sses important for 
galaxy evolution 
(log7VHi^l7.0)ai 



Tripp et al.| ( |2011[ ) show that a LLS 
z = U.927 traces a super-solar metal- 
licity post-starburst galactic wind that extends at least 
to > 68 kpc from its host. On the other end of the 
metallicity spectrum, [Ribaudo et al.| I 
the discovery of a z = Tr274 LLS with 



(]2011b| reported 
he same A^hi as 

in the Tripp et al] study, but with a 2% solar metallicity. 
They show that this absorber resides at impact param- 
eter jO = 37 kpc from a ~ 0.3L* galaxy whose metallic- 
ity is a factor 30 higher than the LLS metallicity. The 
metallicity of this LLS is so low at this redshift that the 
majority of its gas is unlikely to have been recently incor- 
porated into a galaxy. Despite similar A'hi values, the 
metal content, amount of highly ionized gas, kinematics, 
galactic environments, and origins, are all very different 
between these two absorbers. The rest-frame equivalent 
width of the strong Mg II A2796 is VFa = 0.1 and 1 A 
for the metal-poor and metal-rich LLS, respectively; the 
low Mg II equivalent absorber would have been missed 
in most Mg II surveys. 

Prior to these results, only a handful of H I-selected 
LLS or super-LLS (SLLS, 19 < logA^ni < 20.3) 
with inform ation on their m etallicit y were reported at 
z<, 1 (e.g., i Chen fc Prochaska..2000 
Prochaska ct ai.||2(JU4| [Jenkins et al 



Zonak et al .]|2004, 
»; Tripp et al". 



2005| iCooksey et al.||2008| [Lehner et'ar][2009| ) . In [Rib- 



audo et al. l2Ullbl, we argued that the metal-poor LLS 



are the best-yet observational evidence for cold flow ac- 
cretion in galaxies, while the higher metallicity absorbers 
are tracing a mixture of physical processes, such as out- 
flows, recycled gas, and galaxy interactions. Knowing 
the metallicity of these absorbers is critical for placing 
constraints on their origins, but prior to the installation 
of the Cosmic Origins Spectrograph (COS) on the Rub- 
ble Space Telescope (HST), the sample was too small to 
reliably determine the metallicity distribution and the 
range of possible metallicity values in the CGM at z < 1. 
Thanks to the installation of COS, the situation has dra- 
matically improved in just a few years. 

In this paper, we present the first results of an anal- 
ysis of the metallicity of the CGM at z < 1 using LLS 
selected solely based on their H I content. The new LLS 
described in this work were discovered blindly in our Cy- 
cle 17 and in GTO COS programs (see below). In these 



LLS, we can measure accurately the column densities of 
H I and metal ions (in particular, the low ions such as 
C II, Si II, Mg II) that allow us to model the ionization 
conditions and derive the metallicity followi ng a similar 



me thodology used in 
aL] [2008 



survey o: 



>revious studies (e.g., 



Lehner et al 
: H 1-selecte' 



prevK 
12009 



Cookse y et 

Ribaudo et al.[[2011b|. 

ibi. 



provides the first unl 



Uur 

ased 



metallicity study of the LLS (and thus the cool CGM) 
metallicity distribution at z < 1. 

We compare our analysis of the H I-selected with 
metallicity estimates from the literature for SLLS and 
damped Lya absorbers (DLAs, A'hi > 20.3) over the 
same redshift range. However, we emphasize that SLLS 
and DLAs were mostly selected based on knowing a pri- 
ori that they show strong Mg II absorption. Our sample 
of low-z metallicities consists of 28 H I-selected LLS, 29 
SLLS, and 26 DLAs. Seventeen LLS have never been 
published previously, so we have increased the sample of 
LLS by a factor 2 and nearly a factor 5 compared to the 
pre-COS era. 

Our paper is organized as follows. In [2] we describe 
the new and existing LLS samples as wellas the deter- 
mination of t heir metallicities based on photoionization 



models (^2.3 1 with additional technical details (includ- 
ing the description of each absorber) provided in the 
appendix for interested readers. Our main results are 
presented in f|3j where in particular in j |3.1[ we describe 
the bimodal distribution function of the LLS metallicity. 
In this section, we also compare the results from the LLS 
and higher H I column density absorber samples. In Sj4] 
we discuss implications of our new observational results 
for the physical processes and the properties of the CGM. 
Finally, in fjSjwe summarize our main results. 

2. LLS SAMPLE 
2.1. Sample Description 

We search in the UV spectra of QSOs for absorbers at 
z < 1 based on their H I content, i.e., they all satisfy 
the condition 16.2 < logTYni 18.5 within the errors. 
Although the lower TVhi cut-off is somewhat arbitrary, 
it was selected because a break in the continuum at the 
Lyman limit could still be seen in high-signal to noise 
data. These absorbers were found blindly in the UV QSO 
spectra that were obtained for other scientific reasons, 
i.e., in particular, the absorbers were not pre-selected 
based on knowledge of their metal content. But for that 
reason, most of the absorbers have 16.2 < logiVni < 
17.3 where tll < Ij i-G-, the UV flux is only partially 
attenuated. Prior to COS, QSO sightlines were chosen 
to avoid Tll > 1 LLS and very few QSOs were at high 
enough z to produce a large number of LLS. 

We searched for the LLS in four high resolution COS 
programs: 11520 ("COS-GTO: QSO Absorbers, Galax- 
ies and Large-scale Structures in the Local Universe", 
PI Green), 11598 ("How Galaxies Acquire their Gas: A 
Map of Multiphase Accretion and Feedback in Gaseous 
Galaxy Halos", a.k.a. the "COS-Halos" program, PI: 
Tumhnson), 11692 ("The LMC as a QSO Absorption 
Line System", PI: Howk), 11741 ("Probing Warm-Hot 
Intergalactic Gas at 0.5 < z < 1.3 with a Blind Survey 
for O VI, Ne VIII, Mg X, and Si XII Absorption Systems, 
PI: Tripp). Note that targeted systems in the COS-Halos 
program are not in our sample (see a description of these 
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absorbers in 'Tumlinson et al. 2011a Werk et al. 20131 
because the A'hi tor the LLH is not well constrained in 
that sample. At this time, we have not searched system- 
atically the entire COS G130M/G160M QSO archive, 
although this will be underway soon. The redshifts of 
the 12 QSOs are in the range 0.33 < zqso £ 1.44. 
None of the LLS described here are proximate (i.e., with 
Zabs ~ ZQso) systems. 

The choice of the present data sample is motivated 
by three reasons: 1) a large portion of these data was 



already reduced by our team to undertake the original 
scientific goals of these COS programs, which in particu- 
lar helped with the line identification process; 2) several 
of these QSOs were observed with Keck HIRES, allowing 
us to measure Mg II AA2796, 2803 absorption, which is 
a key ion for determining the metallicity and to deter- 
mine the velocity component structure of the absorbers; 
3) most of the spectra have relatively high signal-to- noise 
(S/N> 15-30 per resolution element). 



Table 1 



Name 



log A^Hi 



[X/H] 



IorU 



P 

(kpch-i) 



References 



LLS - 16.2 < logAfni < 19 



J0943+0531 





3542 


16.11 it 0.09 


< -1.30 




< —3.2 


95 


Tholl" 


J 14iy-t-4zU ( 





4256 


iD.i / ± U.Uo 


-1.40 ±0.20 


—2.9 




This paper 


J1435-I-3604 





3878 


16.18 ±0.05 


< -1.40 




> -3.0 




This paper 


PG1116+215 





1385 


16.20 ±0.05 


-0.50 : 




-2.5 


127 


Sem04,Tri98 


PG1522+101 





5185 


16.22 ±0.02 


-0.40 ±0.05 


-3.6 




This paper 


SBS1122+594 





5574 


16.24 ± 0.03 


-1.05 ±0.05 


-3.1 




This paper 


HE0439-5254 





6153 


16.28 ±0.04 


-0.30 ±0.05 


-2.7 




This paper 


PG1338+416 





3488 


16.30 ±0.13 


-0.65 ±0.15 


-3.2 




This paper 


J1419+4207 





5346 


16.34 


-0.20 ±0.20 


-3.9 




This paper 


PG1407+265 





6828 


16.38 ± 0.02 


-1.80 ±0.30 


-2.0 




This paper 


PG1216+069 





2823 


16.40 ± 0.05 


< -1.65 




< -2.6 


3238 


This paper, Prol la 


J1419+4207 





2889 


16.40 ± 0.07 


-0.65 ±0.10 


-3.1 




This paper 


PG1338+416 





6865 


16.45 ± 0.05 


+0.10 ±0.10 


-3.9 




This paper 


PKS0405-123 





1672 


16.45 ± 0.05 


+0.10 : 




-3.1 


100 


CheOO,SavlO 


J1619-I-3342 





2694 


16.48 ± 0.05 


-1.60 + 0.10 


-2.9 




This paper 


PKS0637-752 





4685 


16.48 ± 0.04 


-0.50 + 0.1 


-3.8 




This paper 




n 


2261 


16 61 +012 


-0.80 + 0.25 


-2.8 




Savll 


Ti yi Q 1 Qfin/1 
J i4oo-t-oDU4 


U 


Q'7Qn 

o / oU 


1 a fit; _L n n'z 
ID. DO it U.U ( 


-1.85 + 0.10 


-3.5 




This paper 







7292 


1d.6d ± U.Oo 


< -2.00 




> -3.2 




This paper 


OUT T 





7392 


Id. (I ± U.Oo 


-1.45 + 0.05 


-2.9 




This paper 


r^l\-DU0Oz-D4U 


U 


Q /I c; 1 


1 a nn -L n no 
ID.yU it U.Uo 


< -1.50 




> -3.8 




This paper 




U 




ID.yo It U.UO 


-1.71 + 0.06 


-2.8 


o / 


IXIUIL D 


PGlzUD+4oy 





9270 


T nn _L n 1 n 
1 / .UU ± U.IU 


+0.30 : 




< -2.5 


68 


Trill 


TDT^ci Qno 1 no 


U 


UyoO 


1 7 nn -L n on 
1 ( .UU it U.zU 


< -1.60 




-2.9 


Do 


L^ooUo 


r^Oi0o4+ (UO 


i 


ri/i nn 
U4UU 


1 7 Qn _L n Qn 
1 / .oU it U.oU 


-1.50 + 0.20 


-5.2 




ZjonU4 


DUT 1 Q1 1 


U 


riQ 1 n 


1 7 ciQ -L n nc; 
1 / .yo it U.UO 


-0.19 + 0.08 


-4.0 


Q /I 

o4 


Jenkuo 


PKbUolz- / ( 





2026 


T o oo _L n on 
lo.zz it U.2U 


-0.60 + 0.15 


-3.2 


38 


Leh09 


TON153 





6610 


18.30 ± 0.30 


-1.69 + 0.37 


-3.2 


104 


Kacl2'' 


J1009-I-0713 





3588 


18.40 ± 0.20 


-0.40 + 0.20 


-3.5 


43 


Tumllb,Werl2b 








SLLS 


- 19 < log 


Nni < 20.3 






Q0826-2230 





9110 


19.04 ±0.04 


+0.68 + 0.08 


< -4.1 




Mei09 


Q0005-I-0524 





8514 


19.08 ±0.04 


< -0.47 








Mei09 


Q2352-0028 





8730 


19.18 ±0.09 


< -0.14 








Mei09 


Q1054-0020 





9513 


19.28 ±0.02 


< -0.66 








Des09 


Q2149+212 


1 


0023 


19.30 ±0.05 


< +0.20 








NesOS 


JlOOl-l-5944 





1140 


19.32 ±0.10 


-0.37 + 0.10 


-3.4 




Bat 12 


PG1216+069 





0063 


19.32 ±0.03 


-1.60 + 0.10 


-3.1 


82 


Tri05 


J0928-I-6025 





1538 


19.35 ±0.15 


+0.31 + 0.17 


-3.7 


38 


Bat 12 


J0021-I-0043 





9424 


19.38 ±0.13 


< -0.41 








Des09 


Q1330-2056 





8514 


19.40 ±0.02 


< -0.46 








Des09 


Q1228+1018 





9376 


19.41 ± 0.02 


< -0.37 






38 


Mei09 


Q1009-0026 





8866 


19.48 ±0.05 


+0.24 + 0.15 


-3.7 


35 


Des09 


J1553-I-3548 





0830 


19.55 ±0.15 


-1.10 + 0.16 


-3.2 




Batl2 


J0925-I-4004 





2477 


19.55 ±0.15 


-0.29 + 0.17 


-3.1 


85 


Bat 12 


Q163H-1156 





9004 


19.70 ±0.04 


< -0.15 








Mei09 


Q0153-f0009 





7714 


19.70 ±0.09 


< -0.34 








Mei09 


Q2335-(-1501 





6798 


19.70 ±0.30 


+0.07 + 0.34 : 






Mei09 


J1435-I-3604 





2027 


19.80 ±0.10 


-0.41 + 0.16 


-3.3 


38 


Batl2 


Q2352-0028 


1 


0318 


19.81 ± 0.13 


< -0.51 








Mei09,Des09 


Q0138-0005 





7821 


19.81 ± 0.09 


+0.28 + 0.16 






Per08,Mei09 


J0134+0051 





8420 


19.93 ±0.12 


< -0.36 








Per06 


J1028-0100 





6321 


19.95 ±0.07 


< -0.20 








Des09 


J1028-0100 





7089 


20.04 ± 0.06 


< -0.18 








Des09 


Q1436-0051 





7377 


20.08 ±0.11 


-0.05 + 0.12 




< 33 


Mei09 


Q1455-0045 


1 


0929 


20.08 ±0.06 


< -0.80 








Mei09 


Q1009-0026 





8426 


20.20 ±0.06 


< -0.98 






> 15 


Mei09 


Q1220-0040 





9746 


20.20 ±0.07 


< -1.14 








Mei09 



4 



Lehner et al. 

Table 1 — Continued 



Name 


■^abs 


log A^Hi 


[X/H] logU 


P 


References 






[cm-2] 




(kpch-i) 




J1323-0021 


0.7160 


20.21 ± 0.20 


+0.61 it 0.20 


9 


Per06,Kha04 


Q1107+0003 


0.9542 


20.26 ±0.14 


< -0.51 




Nes08 


DLA - logAfm > 20.3 


TOQOQ 1 nnoo 
J ZoZo-\-v\jZz 


U.DOZU 


on Qo _L n r\'7 
zU.oZ it U.U ( 


— U.4y it U.zz 


■ ■ ■ 


"Di-if nft 
Jr^erUo 




1.0090 


on oii _L nil 
zU.ob ± U.li 


—0. /o ± U.lz 


25 


r^eLUU,ijerSy / ,r^erUD 


C^lbzZ+Zoo 


U.dodO 


on Q£i _L n 1 n 


— 0.8 / ± U.ZO 




OnuUU,ri,aoUU 


C^Iizz-loo 


U.DOzO 


on /I c; _L n n c; 


— 1.00 ± U.lo 




VarUU 


C^lozo-UUzi 


0.7156 


on /I _L n 1 c; 


+0.12 ± U.ZO 




INesUo 


J iDly+oo4i 


O.Oyoo 


on c; _L n 1 n 
ZU.OO ± U.IU 


— U.6o ± U.lo 




Bat 12 


J iDiDH-4i04 


U.ozio 


on fin J- n on 
zU.oU it U.zU 


— U.oO it U.zo 




"D 1 O 

rSatiz 


CjUioy-UUzo 


0.6825 


on art _L n 1 o 
zU.dU ± U.lz 


< —0.08 




INesUo 


Ti nno-Lnvi ^ 




on fifi _i_ n T n 
ZU.Do it U.IU 


— U.UO It U.lo ■ ■ ■ 






/^n A A 1 non 
C^U4o4+Uoy 


O.OOUO 


on (ic\ _L n nil 
zU.oy ± U.Uo 


— U. r y ± U.lz ■ ■ ■ 




PetOO 


J0256+0110 


0.7250 


20.70 it 0.15 


-0.11 it 0.04 




Per06 


Q0256+0110 


0.7252 


20.70 it 0.22 


-0.13 it 0.24 




Nes08 


Q17334-5533 


0.9981 


20.70 it 0.04 


-0.44 it 0.07 




Nes08 


Q2149+212 


0.9111 


20.70 it 0.10 


< -0.93 




Nes08 


Q1229-021 


0.3950 


20.75 it 0.07 


-0.47 it 0.15 


7 


Boi98.Leb97,Ste94 


Q0253+0107 


0.6316 


20.78 it 0.12 


< -0.35 




Nes08 


Q0449-1645 


1.0072 


20.98 it 0.07 


-0.96 it 0.08 




Per08 


J1107+0048 


0.7400 


21.00 it 0.03 


-0.54 it 0.20 




Per06,Kha04 


Q1007+0042 


1.0373 


21.15 it 0.24 


-0.51 it 0.24 




Nes08 


Q1727+5302 


0.9445 


21.16 it 0.05 


-0.56 it 0.06 




Nes08 


J1431+3952 


0.6019 


21.20 it 0.10 


-0.80 it 0.20 




E1112 


Q1328+307 


0.6920 


21.25 it 0.06 


-1.20 it 0.09 


6 


Mey95,Leb97,Ste94 


Q1225+0035 


0.7728 


21.38 it 0.12 


-0.78 it 0.14 




Nes08 


Q1727+5302 


1.0306 


21.41 it 0.03 


-1.36 it 0.08 




Nes08 


2353-0028 


0.6043 


21.54 it 0.15 


-0.92 it 0.32 




Ncs08 


Q0235+164 


0.5240 


21.65 it 0.15 


-0.58 it 0.15 


6 
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i98 , Bur96: Burbidgc ct al. 1996| (JheO O: |(Jhen fc P Fochaska 
2000 OhuOO: Churchill ct al. 2000 T^oOS: IC ooksey et 'al. 2008 
es09: Dcssaugcs-Zavadsky ct al. 2009 ElUz: Ellison ct al. 2012 
Kacl2 : Kacprzak ct al. 20i2b Jcn05: Jenkins ct al. 2005, Kha04: 
Lc Brun ct al. 1997 Lch09 : |Lehner | 



Khare et al. 



Lcb97: 



et al.||ii(iOO'rMei09: Mciring ct al. 2009 Mcy95: Meyer et al. 
raW^NesOS: 'Nestor ct al. 2008 Fer06: Pcroux ct al. 200(7 Per08:' 
Pcroux ct al. 2008, PctOO: Pcttini ct al. 2000 Prolla: Prochaska 
et al. 2011a RaoOO: Rao & Turnshek 2000, Riblla: Ribaudo et 
al. 2011b SavlO: Savage ct al. 2010 Savll: Savage ct al. 20f||' 
^em04: .Sembach ct al. 2004 Stc94: Stcidcl et al. 1994, TholE 
I'hom et al.||2Ull| Tri98: 



iLlli^OllI T ri98: Tripp et al. 1998 Tri05: Tripp ct al 

ii005| 'iHlliTi'ripp et al 2011 I'uml 1 b: .I'umlinson et al- , .^01 1 b 



[ru rnshek ct al.|2l)04| !^on04: [Zonak el al.p0(j4| V('erlijb: 
IWerk et al. 2013 

Note. — We express the abundances with the usual logarithmic 
notation. U is the ionization parameter derived from the Cloudy 
simulations and p is the impact parameter to generally the closest 
gala xy (there might be evidence of other galax ies within 150 kpc, 
e.g., [Jenkins et al.|2005l iKacprzak et al.|2012b|| For the LLS, the 
metallicity is estirnated irom a elements (e.g, U, Mg, Si). For the 
SLLS and DLA, the metallicities were estimated from Zn, except 
those estimated by Batl2 and Tri05, where X is an a element. For 
the SLLS with logA'^jjl $ 19.7 and no U value, no ionization cor- 
rection was estimated. For these absorbers, ionization corrections 
may be necessary, and we emphasize this uncertainty by adding 
a colon in the metallicity column and adopt the me asured value 
as an upper limit, a: this absorber was analyzed by [Thorn et al.| 
|2011| but no error was given on A^h i and la upper limits were 
used to estimate the metallicity. Using H I AA926, 920, 919, 918, 
we re-estimated A'^hi; we also re-estimated the metallicity using 
3(7 upper limits. Our and Thom et al. results are consistent within 
Ict. b: As this absorber was pre-selected on Mg II absorption it is 
not included in our sample of H I-selected LLS. 

To complement our sample of unpublished LLS, we sys- 
tematically surveyed the literature for H I-selected LLS 
at z < 1 for which the metallicity has been estimated. 
Our entire LLS sample is summarized in Table [T] where 
the redshift of the absorbers, the H I column density, the 



metallicity, ionization parameter (see below), and the im- 
pact parameter (p) are listed. We assume the absorber is 
associated with the galaxy near the QSO with a redshift 
consistent with that of the absorber. In some cas es, mul- 
tiple galaxies may be wit hin abo ut 100 kpc (e.g., Jenkin s 



_ga 
11 



et al. 2005 



In cases w 



Cooksey et al. 2008 ; Tumlinson et al 
lere more than one galaxy is observed 



the 



galaxy with the smallest impact parameter p and near- 
est redshift was selected. We emphasize that the current 
galaxy redshift surveys are quite heterogeneous in terms 
of depth and completeness, and it remains possible, e.g., 
that there are very low luminosity galaxies closer to the 
sightline or associated luminous galaxies. 

To the best of our knowledge all LLS but one summa- 
rized in Table [l] were discovered serendipitously in the 
HST UV spectra, i.e., they were not pre-selected owing 
to existing optical observations of Mg II absorption. For 
most of the QSOs observed prior the COS era, a member 
of our group was involved in the original scientific pro- 
gram that led to the acquisition of the QSO. For the other 
two QSOs (PG1634+704 and TON153), we searched the 
MAST archive for the abstract of the original program to 
make sure that the program was not pre-selected based 
on the strength of a metal line. The only ex ception is the 
metal- poor LLS observed toward TON153 (Kacprzak et 
aL][2012b| , which was originally targeted for its known 
weak {W ~ 0.3 A) Mg II absorber at z = 0.6610. The 
number of H I-selected LLS is 28, while our total sample 
of LLS at z < 1 consists of 29 LLS. 



2.2. New LLS Sample 
2.2.1. Data Reduction and Presentation 
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The 17 new absorbers were all observed with COS. In- 
formatio n on the design and performance of COS can be 



found in Green et al. 



(2012). The data reduction and 



co-addition p r ocedu res ot the COS data are de scribed in 



Thorn et al.| ( |2011[ ) and [Meiring et al.| ( [2011 ), and we 

refer the reader to these papers tor more details on COS 
data processing. In short, the individual exposures were 
co-added in photon counts and the shifts between the 
individual exposures were determined using the Galactic 
interstellar absorption lines or well-detected lines from 
multiple-line extragalactic systems. In either case, align- 
ment was based on comparison of transitions of compara- 
ble strength. The exposures were shifted to the common 
heliocentric frame and co-added. The resolution of the 
COS G130M and G160M gratings is i? « 17,000. The 
redshift was determined from the peak optical depth of 
the strongest H I component using the weak H I transi- 
tions. When they are detected, we find that C II, O II, 
and Mg II are at the same redshifts as the strongest com- 
ponent of H I. We normalized the QSO continua with low 
order Legendre polynomials near the absorption under 
consideration. Typically polynomials with orders m < 3 

were fitted. 

In the Figs. |A2| to |A18| of the Appendix, we show for 
each LLS the normalized profiles of H I (generally a weak 
and a strong transition) and all or most of the observed 
metal lines covered by COS as a function of the rest- 
frame velocity. We also show the normalized profiles of 
Mg II when available. For all the QSOs, except for two 
sightlines near the LMC (PKS0552-640, PKS0637-752) 
and the two GTO QSOs (HE0439-5254, SBS1122-f594), 
high-resolution (about 6 knis^^ FWHM) optical spec- 
tra were obtained with HIRES on the Keck I telescope. 
For the LMC sightlines (program 11692), optical spectra 
were obtained with the MagE spectrograph on the Mag- 
ellan Clay telescope at Las Campanas Observatory. The 
MagE resolution is much cruder than Keck HIRES res- 
olution, about 55 kms~^, but sufficient to constrain the 
total equivalent width and column density of Mg II. In- 
formation about the reduction and normali zation of the 
Keck and Magellan da t a can be found in |Thom et al.| 



( |2011[ ) and |Werk et al.| ( |2012[ ), respectively, 
because the Mg II AXFTW, 2803 doublet is so strong, 
we will demonstrate below that it is an excellent metal- 
licity tracer for the LLS when coupled with A'h i derived 
from the UV data and within the studied A'hi inter- 
val. The high resolution Mg II spectra also show that 
its component structure is relatively simple for most of 
the absorbers in our new LLS sample. In general one 
component largely dominates the column density pro- 
file. The exception in the new sample is the absorber 
z = 0.2889 toward J1419+4207 where there are at least 
3 components, two of them dominating the column den- 
sity profile. 

2.2.2. Column Densities of the Metal Ions 



Table 2 

Adopted Column Densities" 



Species 


logAf 
[cm~^] 


■"1, f2 
(kms~^) 





Table 2 — Continued 



Species log A'' vi,V2 

[cm^^] (kms^^) 



PHL1377 - 2abs = 0.7392'' 

lil 16.72 ±0.03 -82,100 

CII 13.35 ±0.08 -42,45 

cm > 14.11 -82,100 

CIV 14.17 ±0.05 -82,100 

Oil < 14.30 ±0.04 -82,100 

OIII > 14.83 ±0.02 -82,100 

OIV > 14.90 ±0.03 -82,100 

OVI 14.22 ±0.11 -82,100 

Mgll 11.98 ±0.07 -42,45 

Si IV 13.25 ±0.09 -82,100 

SIII 13.17 ±0.05 -82,100 

SIV < 13.30 -82,100 

SV 12.68 ±0.10 -82,100 

SVI < 13.06 -82,100 

AlII < 12.10 -42,45 



PG1338+416 - Zabs = 0.3488 

in 16.30 ±0.13 -40, 75 

CII 13.90 ±0.05 -40,75 

cm > 14.12 ±0.07 -90,75 

OVI 14.56 ±0.04 -90,75 

Mgll 12.49 ±0.02 -40,75 

Sill 12.65 ±0.15 -40,75 

Sim > 13.21 ±0.06 -40,75 



PG1338+416 - Zabs = 0.6865 

in 16.41 ± 0.02 -75,102 

CII 14.21 ±0.05 -45,102 

cm > 14.11 -57,125 

Oil 14.30 ±0.04 -45,102 

OIII > 14.92 -75,125 

OIV > 15.10 -75,125 

OVI 14.75 ±0.05 -75,135 

Mgll 13.12 ±0.04 -35,75 

SIII 13.52 ±0.08 -35,75 

SiV 13.58 ±0.04 -75,110 

SV 12.94 ±0.06 -75,110 

SVI 13.22 ±0.06 0,100 



PG1407+265 - ^abs = 0.6828 

in 16.38 ±0.02 -75,80 

CII < 12.75 -75,80 

cm > 13.97 ±0.04 -75,80 

Oil < 13.18 -75,80 

OIII > 14.65 ±0.01 -75,80 

OIV > 14.87 ±0.03 -75,80 

OVI 14.07 ±0.03 -75,80 

Mgll < 11.47 -75,80 

SIII 13.11 ±0.04 -75,80 

SiV 13.17^0 24 -75,80 

SV 12.80 ±0.05 -75,80 

SVI 12.73 i°i3 -75,80 



PG1522+101 - Zabs = 0.5185 

in 16.22 ±0.02 -60, 50 

CII 13.37 ±0.05 -20,50 

cm 13.23 ±0.04 -30,70 

Oil 14.08 ±0.04 -30,30 

OIII 13.73 ±0.06 -20,60 

OIV < 13.49 -30,30 

O VI < 13.34 -30, 30 

Mgll 12.28 ±0.01 -25,25 



PG1522+101 - Zabs = 0.7289 

III 16.66 ±0.05 -60,60 

CII < 12.80 -60,60 

cm 13.40 ±0.15 -30,60 

Oil < 13.42 -60,60 

OIII 13.88 ±0.05 -50,60 

OVI 13.92 ±0.07 -30,40 
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Table 2 — Continued 



Species 


log AT 
[cm-2] 


Vl,V2 

^Km s ) 


Mg II 
S III 
S IV 
S V 


< 11.49 

< 12.44 

< 12.70 

< 12.28 


-60,60 
-60,60 
-60.60 
-60,60 


PKS0552-640 - Zabs = 




H I 

C III 
O VI 
Mg II 
Si II 
Si III 


16.90 ±0.08 
12.92 ±0.07 

< 13.24 

< 12.22 

< 12.27 
12.41 ±0.09 


-80.80 
-30^60 
-60,60 
-60.60 
-60,60 
-60,60 


PKS0637-752 - z^bs = 


0.4685 


H I 
C II 
C III 
O II 
O III 
O IV 
O VI 
Mg II 

s IV 
s V 


16.48 ±0.04 
13.72 ±0.03 

> 13.82 ± 0.03 
14.03 ±0.04 

> 14.45 ± 0.04 : 
14.32 ±0.05 
13.87 ± 0.05 
12.82 ±0.04 

< 13.40 
12.61 ± 0.08 


-60,60 
-50,50 
-90,60 
-50,50 
-90,60 
-90.60 
-80.60 
-60,60 
-90.60 
-90,60 


HE0439-5254 - Zabs = 


0.6153 


H I 
C II 
C III 
O II 
O III 

o IV 
o VI 

S III 
S IV 

s V 
s VI 


16.28 ±0.04 
13.93 ±0.03 

> 14.49 
14.19 ±0.03 

> 14.95 : 

> 15.27 : 
14.84 ± 0.03 

13.68 ±0.04 

13.69 ±0.06 
13.24 ± 0.04 
13.13 ±0.10 


-45,70 
-45,50 
-81,89 
-45,50 
-112,50 
-112, 70 
-150,210 
-45.60 
-45,60 
-45,60 
-45,60 


SBS1122+594 - Zabs = 


= 0.5574 


H I 
C II 
C III 
O II 
O III 

o IV 
o VI 

S IV 


16.24 ± 0.03 

< 13.60 
12.99 ±0.08 
13 61+"-^° 
13.90 ±0.08 

< 13.43 
13.99 ±0.12 

< 13.49 


-60,40 
—40, 40 
-40, 40 
-40,40 
-40,40 
-40,50 
-40,40 
-40,40 


J1419+4207 - Zabs = 


0.2889 


H I 
C II 
C III 

o VI 

Mg II 
Si II 
Si III 


16.40 ±0.07 
13.92 ±0.13 
> 14.00 ±0.08 
14.54 ± 0.04 
12.58 ±0.02 

< 12.96 
13.33 ±0.08 


-80.90 
-80,90 
-80,90 
-80,90 
-30,40 
-80,90 
-80,90 


J1419+4207 - Zabs = 




H I 
C II 
C III 
O VI 
Mg II 
Si III 


16.17 ± 0.06 

< 13.32 
13.33 ±0.10 
13.76 ±0.11 

< 11.78 
12.85 ±0.13 


-70,70 
-50.60 
-50,60 
-80,20 
-50,60 
— 50. 60 


J1419+4207 - 2abs = 


0.5346 


H I 
C II 
C III 
O II 


L0..3t _0.12 

< 13.64 
13.56 ±0.09 
13.97 ± 0.11 


-40,50 
-40, 50 
-45,50 
-40,50 



Table 2 — Continued 



Species 


log AT 






[cm^^] 


(kms^^) 


O III 


< 14.04 


-40, 50 


O VI 


< 13.90 


-40, 50 


Mg II 


12.64 ±0.02 


-34,30 


J1435+3604 - Zabs = 


0.3730 


H I 


16.65 ±0.07 


-70,60 


C III 


13.15 ±0.09 


-70,60 


O II 


< 14.02 


-70,60 


O III 


< 14.31 


-70,60 


O VI 


< 13.51 


-70,60 


Mg II 


11.56 ±0.12 


-20,20 


Si II 


< 12.54 


-70,60 


Si III 


< 12.58 


-70,60 


J1435+3604 - Zabs = 


0.3878 


H I 


16.18 ±0.05 


-70,60 


C II 


< 13.64 


-70,60 


C III 


> 14.01 ±0.07 


-130, 100 


O II 


< 13.99 


-70,60 


O III 


> 14.45 ±0.13 


-110,30 


O VI 


14.24 ±0.09 


-110,30 


Mg II 


< 11.66 


-70,60 


Si II 


< 12.69 


-70,60 


Si III 


< 12.71 


-70,60 


J1619+3342 - Zabs = 


0.2694 


H I 


16.48 ±0.05 


-60, 70 


C II 


< 13.37 


-60,60 


C III 


13.65 ±0.05 


-60,60 


O VI 


< 13.93 


-60,60 


Mg II 


11.75 ±0.12 


-40, 40 


Si II 


< 12.54 


-60,60 


Si III 


13.02 ±0.03 


-60,60 


Si IV 


< 12.90 


-60,60 


PG1216+069 - Zabs = 


0.2823= 


H I 


16.40 ±0.07 


-60, 70 


C II 


< 13.15 


-40, 40 


C III 


> 13.63 


-35, 40 


N II 


< 13.03 


-25,25 


N III 


< 13.03 


-25,25 


O VI 


13.60 ±0.05 


-25,50 


Si II 


< 11.85 


-25,25 


Si III 


> 13.10 ±0.09 


-25,25 







Note. — "The reader should refer to the Appendix for a de- 
scription of each of these absorbers (in particular for the complete 
information about the velocity structure; the velocity interval used 
for the integration of the Na{v) profile for each species is listed in 
the last column. ''C IV, Si IV, and Al II are from STIS E230M 
data. <=The analysis of PG1216+069 is based on both COS and 
STIS E140M data; see Appendix for details. The symbol ">" in- 
dicates that the absorption profile is saturated. The symbol ">" 
indicates that the absorption profile might be saturated. The sym- 
bol "<" indicates that there is no detection at the 3cr level. The 
symbol "<" indicates that the profile might be contaminated by 
an unrelated absorber. A value followed by ":" indicates that it 
is uncertain (e.g., a line is saturated, but could also be contami- 
nated) . 

We employed the apparent optica l depth (AOD) method 
described by 'Savage & Sembach ( 1991 ) to estimate the 
column density of the metal ions. I'he absorption profiles 
are converted into apparent column densities per unit 
velocity, Na{v) = 3.768 x WHn[F,{v)/FoUv)]/{fX) 
cm^^ (kms^^)^^, where Fc{v) and Fobs(w) are the 
modeled continuum and observed fluxes as a function of 
velocity, respectively, / is the oscillator strength of the 
transition and A is the wavelength in A. The atomic pa- 
rameters are for the UV transitions from Morton ( 2003 1 
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and for the EUV transitions from [Verner et a~ (1994). 
The total column density was obtained by integrating 
over the absorption profile Na — J^^ Na{v)dv. When no 
detection was observed, we estim ated a 3g upper lim it 
followi ng th e meth od described by [Lehner et al. ( 2008 ) . 

Figs. |A2| to [AT8| in the Appendix show that often sev- 
eral transitions of the same ion are available to estimate 
the column density (e.g., for O II, O III, O VI, S III, 
respectively). This allows us to assess if the lines are 
saturated and/or contaminated. For the Mg II doublet, 
the only evidence of saturation is for the z = 0.6865 ab- 
sorber toward PG1338-f416, but even in this case the 
saturation is mild since the difference in the Na values 
between the strong and weak transitions is small (imply- 
ing a correc tion of only 0.10 dex follo wing the method 
described in iSavage & Sembach 1991). Other doublets 
(e.g., O VI) or lines with multiple transitions (e.g., O II) 
generally give similar integrated column densities within 
Icr error, except when there is evidence for contamina- 
tion. For assessing the line saturation in species with 
a single transition (e.g.. Si III, C III), we use the infor- 
mation from the profiles of other species with multiple 
transitions that trace similar type of gas (i.e., where the 
temperature and turbulence should be similar). In Ta- 
ble [2] we summarize our adopted column densities and 
the velocity interval for the integration of the profi les. 
W hen m ore than one transition is available (see Figs. A2 



to |A18 ), the adopted column density represents the av- 
erage value between the different estimates when there is 
no evidence of saturation or contamination (i.e., the Na 
values agree within Icr error and there is not a systematic 
increase of the Na with decreasing /A). 

2.2.3. HI Column Densities 

We estimate the H I column density from the weak 
transitions of the Lyman series down to A916 (typically, 
the Lyman series lines become too blended and/or the 
continuum placement too uncertain at A < 916 A to es- 
timate reliable A^hi or W\) and from the flux decrement 
at the Lyman limit. In the latter method, by estimating 
the optical depth tll at the Lyman limit (All ~ 912 
A), we can directly deduce A^hi since A'hi — tll/cll 
where ctll = 6.3 x 10~^^ cm^(A/912)'^ is the absorption 
cross section of the hydrogen atom at the Lyman limit 
( Spitzer|1978 ). To define the continua of the QSO before 
an d after the break, we a dopted the method described 



HE0439-5254 (Fig. AlO), Ly/3 has several components 
at u > 75 kms""'^ and v < —75 kms~^. The strongest 
H I component for which log A'hi ^ 16.2 is observed over 
— 75 < < 75 kms~^, and its column can be estimated 
using only the weaker transitions (e.g., H I AA923, 919). 
This ensures the metallicity estimates are accurate, since 
the low ions are generally seen in the strongest H I com- 
ponent only. 

We also use a third method where we measure the 
equivalent width {W\) for all the Lyman transitions that 
do not contain any evidence of widely displaced compo- 
nents to estimate A^hi from a single-component curve- 
of-growth (COG). For example, for the absorber toward 
PHL1377, we use both strong and weak lines for the 
COG, while for the absorber toward HE0439-5254, we 
only use the weak transitions as there are several com- 
ponents in the stronger ones both blueshifted and red- 
shifted with respect to the main principal components. 
The COG method used a mini mization and Y^_er 



(1990). 



y m es- 



ror derivation approach outlined by Savage et al. 
The program solves for log N and b independent 
timating the errors. 

In the Appendix, we summarize for each absorber 
which methods were used to determine A'hi. In gen- 
eral, there is a good agreement between these methods 
when we were able to use two or three of them. The 
three methods explore different parameters and contin- 
uum models, allowing us to reliably estimate A'hi and 
its error. When more than one method was used, we av- 
eraged the values and propagated the errors accordingly. 

2.3. Determination of the Metallicity of the LLS 

Large ionization corrections are required to determine 
the metallicity from the comparison of H I with singly 
or doubly ionized metal species, except where O I is 
observed, which only happened for the three strongest 
LLS in the s ample (with log Ah i > 18). As displayed 



in Figs. A2 to A18 (see also references given in Tablejl]), 



by Ribaudo et al. (1201 fal) by fitting a composite QSO 



the ionization structure can be quite complicated in some 
cases, with the detection of species in different ionizing 
stages and sometimes with velocity shifts between the 
different ions. However, in other cases, it is much simpler 
with very few ions detected, and, as we will see below, 
these generally correspond to low metallicity absorbers. 
While it is evident that several ionization processes are 
at play for some of the absorbers (see below), our main 
aim here is not to understand the detailed origins of the 
physical processes that give rise to the observed ioniza- 
tion and velocity structures (we will explore this further 
in the future). Instead our goal is to determine the metal- 
licity of the weakly ionized gas that is observed over the 
same velocity interval as the strongest H I component 
with logA^Hi ^ 16.2. This requires we concentrate on 
modeling the ionization of the singly and doubly ionized 
species. 



spectrum developed by Zheng et al.| ( |1997[ ) . Depending 
on the complexity of the QSO continua and interstellar 
spectra (in particular near the Milky Way Lya absorp- 
tion), this method could be reliably used to determine 
A'hi for absorbers with log A^hi ^ 16.3 or tll ^ 0.15 for 
the typical S/N levels in our sample. 

To estimate the H I column density from the Lyman 
series, we employed the AOD method described above. 
We use all the available weak Lyman transitions that 
are ot contaminated and for which the continuum can 
be reliably modeled. As for the metal lines, we use 
several H I transitions to assess the saturation level, 
if any. As we are solely interested in the components 

with log Ah I > 16.2, we only integrate the Na{v) pro- scribed by [Ferland et al.|[T998[ ). We assume the gas 
files over the velocity interval that defines these compo- 
nents. For example, for the z = 0.6153 absorber toward 



2.3.1. Photoionization Models: Methodology 

To determine the metallicity and physical conditions 
of the strong H I component in each absorber, we 
model its i onization using Cloud y (version clO, last de- 



is photoiomzed, modeling it as a uniform slab in ther 
mal and ionization equilibrium. We illuminate the slab 
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Figure 1. Distributions of A^hi (left) and the Mg II A2796 rest-frame equivalent width (right) for the metal-poor and more metal-rich 
absorbers. Only absorbers from the H I-selected LLS sample are considered. The hashed histograms highlight values that are upper limits, 
i.e., the strong transition of the Mg II doublet is not detected. Visually and statistically, the differences between the metal-poor and more 
metal-rich samples are not significantly different for the A'^hi distribution, but they are for the W\{TsAg II) distribution. 
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Figure 2. The column densities of Mg II and O VI versus the 
metallicity of the cool, photoionized gas for the LLS in our sample. 
The A^o VI versus [X/H] is essentially a scatter diagram, but there is 
a strong correlation between A'^Mg II and the metallicity. The dotted 
line show a linear fit to the data with error bars using a "robust" 
least absolute deviation method (logA^Mgll = 0.9 [X/H] -|- 13.2). 
The departing point at low metallicity is the Mg Il-selected LLS 
with a large A'jjj (18.3 dex) compared to the other LLS (16.1-17.1 
dex). 



with the Haardt-Madau background radiation field from 
quasars and galaxies (H1VI05, as implemented within 
Cloudy). We describe in detail two examples of Cloudy 
photoionization models in the Appendix. In summary, 
for each absorber, we vary the ionization parameter, 
U = nj/nn =H ionizing photon density/total hydrogen 
number density (neutral + ionized), and the metallicity, 
[X/H] = log Nx/Nu - log(X/H)Q, to search for models 



that are consistent with the constraints set by the column 
densities determined from the observations. W e assume 
solar r elativ e heavy element abundances from |Asplund] 
et al.| ( 2009 ) . We do not include a priori the effects of 



dust or nucleosynthesis on the relative abundances, al- 
though we consider these possibilities a posteriori. The 
preferred species to constrain the ionization parameter 
are those where their velocity structures follow the H I 
velocity profile the best. The velocity structure of the 
singly ionized species (e.g., O II, C II, Si II, and Mg 11) 
match well the H I velocity structure, and, in some cases, 
doubly ionized species can also be used reliably based on 
their similar velocity profiles compared to the low ions 
and H I. The excellent correspondence between the peak 
optical depths of the H I and the strongest low-ion lines 
implies the H I and low ions must be co-spatial. This 
gives us confidence that the metallicity distribution re- 
sulting from this modeling is not a simple consequence 
of unrelated absorption. However, that is why we do not 
blindly compare the total column densities: various ve- 
locity components may have different metallicities, and 
While the total H i column density may not change by 
more than ^-^0.2-0.3 dex compared to the TVhi value in 
the strongest component, that may not be the case for 
metal ions such C iii and O iii, where most of the column 
density may in fact be assoc iate d wit h the low A'h i com- 
ponent (see, e.g., Figs. A6 and All). When singly and 
doubly ionized species can be used, ionic ratios of the 
same element can be used to determine U and [X/H], 
alleviating uncertainties from possible non-solar relative 
abundances. When singly ionized species are not de- 
tected, the 3tT upper limit on N-^g ii provides a stringent 
constraint on the metallicity (see next), and we also use 
doubly ionized species (C III, O III, Si III) if their velocity 
profiles follow well the profile of the strongest H I com- 
ponent. The metallicities of the LLS were determined 
using a elements (O, Si, Mg, S). 

For the LLS in our sample, low ionization species like 
Mg II are powerful metallicity indicators, as we illus- 
trate in the right panel of Fig. [T] and the top panel of 
Fig. [2] In Fig. [T] we show the distribution of the rest- 
frame equivalent width, W\, of Mg II A2796 color coded 
according to the metallicity of the absorbers: metal-poor 
absorbers with [X/H < — 1] have a relatively weak Mg II 
A2796 absorption (Wx < 100 mA, with Wx < 30 ml 
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for more than half of the weak Mg II absorbers), while 
more metal-rich absorbers have typically strong Mg II 
absorption {W\ > 100 niA, this figure does not include 
the extremely strong Mg II absorber - Wx '--^ 1 A - with 
supersolar metall icity toward PG1206-I-459; see Table [T] 
and [Tripp et aL]|201l[ ). Fig. [2j demonstrates further a 
strong correlation between [X/H] and-/VMgii- This strong 
metallicity dependency is a result both of the limited 
range of TVhi (mostly 16.2 < log A'hi ^ 17) and narrow 
log U interval (see Appendix) sampled by our absorbers. 
The column densities of the doubly ionized species (like 
C III or Si III) can, however, vary by a large factor over 
the constrained U interval. 

The LLS that are already in the literature were mod- 
eled following a similar method in the original publica- 
tions (see Table IT] for the references). For a very few 
absorbers (especially pre-2 005), the HM96 (instea d of 
the HM05) UV background ( |Haardt fc Madau|1996[) was 
adopted, and sometimes the UV background from QSOs 
only was favored. This may introduce so me differences 
in the metallicity (e.g., Howk et al. 20091, but no more 
than a factor about 2-3, and, specincaily for the strong 
H I absorbers, J. Werk et al. 2013 (in preparation) show 
that the impact of changing the UV background is quite 
small on the derived metallicities. Further, as this con- 
cerns only four systems, it would not affect our conclu- 
sions anyway. Because most studies also use principally 
Si and Mg to constrain the metallicities of the photoion- 
ized gas, the change in the solar abundances over the last 
8-10 years has little effect since the changes are minor 
for these elements. 

In the Appendix, we discuss the results of our pho- 
toionization modeling for each absorber (and see the ref- 
erences given in Table [T] for the absorbers that are al- 
ready in the literature). The adopted metallicities are 
summarized in Table IT] The errors on the metallicity 
and U reflect the range of values allowed by the Icr un- 
certainties on the observed column densities. They do 
not include errors from the limitation of the models used 
to estimate the ionization corrections. 

2.3.2. Photoionization Models: Systematic Uncertainties 

For all the LLS in our sa mple , the fraction of H that 
is ionized exceeds 90% (see [3.3 1, and the metallicity es- 



timates are therefore photoionization model dependent. 
We now assess whether we can identify subtle ionization 
effects that may produce spurious metallicity systemat- 
ics and whether trends can be observed between the ob- 
served absorber properties and the metallicity that might 
suggest unidentified systematics. As we discuss in §3.11 
we find that the metallicity distribution of the LLS is 
bimodal, and it is important to demonstrate that this 
result is robust and not an artifact from the large ioniza- 
tion corrections required to estimate the metallicity. 

Since the absorbers straddle the boundary of the opti- 
cally thin and thick regimes, there might be some depen- 
dence of the results with TVhi- In Fig. [l] (left panel), we 
show the A^Hi distribution for two metallicity intervals, 
larger or smaller than 10% solar (see Sjs] for the justifi- 
cation of these two metallicity intervalsjT Both visually 
and statistically, the differences between the metal-poor 
and more metal-rich samples are not significant. There 
is no dependence of our derived metallicities on A^hi (see 



also the Appendix). 

In order to reduce potential systematics from differ- 
ent ionization processes, we (and others, see references 
in Table [T]) typically do not rely on higher ionization 
stages (e.g., comparing O II and O IV) to constrain the 
ionization models. High ions can be produced both by 
photoionization and coUisional ionization processes and 
disentangling these processes can be quite complicated 
and can lead to erroneous metallicity estimates. For ab- 
sorbers with 16.2 < logA^Hi ^ 18.5, the photoionization 
models that match the properties of the low ions system- 
atically underestimate the column densities of the high 
ions (see Appendix and references in Table [T]). The inter- 
pretation of this discrepancy is that the low and high ions 
are not cospatial and that other physical processes are at 
play. The fact that we do not rely on the high ions to 
constrain our models also has the benefit of reducing the 
uncertainties from leaking phot ons from galaxies near the 
absorbers: for the Milky Way, [Fbx et al. (2005) showed 
that C II/C III (Si Il/Si III) for gas at ~~5irkpc from 
the Milky Way is almost negligibly affected if the gas is 
immersed in the extragalactic UV background alone or 
in the UV background plus the Galactic radiation field. 
Larger effects are typically observed for C II/C IV (or 
Si Il/Si IV) between the UV background alone and the 
UV background plus the inclusion of the local radiation 
field because the Galactic field does not include many 
photons above the ionization edge. 

We minimize dust or nucleosynthesis effects by com- 
paring when possible single elements in different ioniza- 
tion stages (e.g., C II/C III, O II/O III). We also reduce 
the nucleosynthesis effects by considering when possible 
ratios of a elements (O, Mg, Si, S), as their relative abun- 
dances should remain close to solar values (except in the 
presence of strong dust depletion, which should be neg- 
ligible for the LLS in our sample as discussed below). 
On the other hand, a large fraction of C is produced 
in intermediate-mass stars, while, e.g., O is thought to 
be synthesized almost entirely in massive stars. This 
leads to a time lag in the production of these elements 
and hence [0/C] = log A^o/A^c - log(O/C)0 may not be 
zero. Instead, in galaxies, [C/0] ~ —0.4 to —0.2 when 
[O/Hl < —0.5 (with a large scatter, e.g., [Garnett et aL 
1999; [Henry et al. 2000, and references therein). In the 
diffuse environment that the LLS probe, dus t depletion 
is expected to be small or absent (but see Menard & 



Fukugita 2012) for the species considered in this work 
(in particular, C, O, Mg, Si, S) and many of the species 
(C, O, Mg, Si) are depleted by a similar amount in "halo" 
condition based on the r esults from gas depletion studies 



in the M ilky Way halo (Savage & Sembach 1996 
ins||2009|) . Furthermore, in the Galactic halo, C, CT 
Si are all depleted by the 



Jenk- 

ML, 

depleted by the same amount. Our Cloudy 
photoionization models do not suggest strong departures 
from non-solar relative abundances for O, Mg, Si, S. Nev- 
ertheless it is possible that there could be a factor ~2 
departure from solar relative abundances, which could 
then affect by a similar factor our metallicity estimates. 
However, while this effect may shift the mean metallicity 
in both branches of the metallicity distribution, the ef- 
fect would be too small to alter the bimodality nature of 
the metalhcity distribu tion described in ^3.1 For [C/a], 
we refer the reader to §3.2| where this ratio is discussed 
in more detail. 
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Figure 3. Metallicity distribution function of the LLS at z < 1. The hashed histograms highlight values that are upper limits. The 
metallicity distribution is bimodal where the peak values are shown by the vertical dotted lines. For the LLS with [X/H] < — 1, the mean 
metallicity is an upper limit as 6/14 of the [X/H] estimates are upper limits. 



We also checked whether coUisional ionization mod- 
els could explain the observatio ns a s well as the pho- 
toionization models discussed in 92.31 We used the coUi- 



si onal ionization eq uilibrium and non-equilibrium models 
of Gnat fc Ster nberg (2007) and shock ionization mod- 
els oT Gnat fc Sternberg (2009|) to compare the predicted 
and observed ionic ratios ot (J II/C III, Si Il/Mg II, and 
O ll/Mg II. For the majority of the absorbers, the ob- 
served ionic ratios are well outside the predicted curves 
by these models. For 2 absorbers, only solar or super- 
solar metallicity shock ionization models may be ade- 
quate (for these, only limits on the ionic ratios are avail- 
able). Hence the low ions are best explained by pho- 
toionization models. 

We therefore conclude that there does not seem to be 
any hidden artifact from the ionization correction s on 
the metallicity distribution of the LLS discussed in §3.1[ 
We did not find any relation between the metallicity and 
A'hi (see also appendix), which spans a factor > 200 and 
160, respectively, between the lowest and highest values. 
On the other hand, the photoionization models naturally 
explain the strength of the observed absorption profiles 
for the singly ionized metal species. 

3. PROPERTIES OF THE LLS 

3.1. Metallicity Distribution Function 

Fig. [3] shows the metallicity distribution function for 
the 28 H I-selected LLS at z < 1 summarized in Table [l] 
Visually the distribution is strikingly bimodal. The up- 
per half of the sample (14/28) has a mean and dispersion 
([X/H]) = -0.33 ± 0.33 (median -0.40). The lower half 
has ([X/H]) < -1.57±0.24 (median < -1.60); the mean 
and median are strictly an upper limit because only an 
upper limit could be derived for the metallicity of 6/14 
absorbers. The above is based on a visual inspection of 



the histogram distributions of the LLS metallicity. Un- 
surprisingly, a K-S test shows that the LLS metallicity 
distribution is not consistent with a normal distribution. 

To further quantify whether the LLS metallicity distri- 
bution is unimodal or bimodal, we use a gaus sian mix- 
ture model i ng (G JMIVI) algorithm developed by Muratov 



& Gnedin (2010) as well as an independent test ot the 
bimodality based on the dip statistic. These statisti- 
cal tests were used to quantify and test the bimodal- 
ity of the metallicity distribution of the globular clusters 
in the IVIilky Way, which also shows two peaks in the 
metallicity. The dip test does not make any assump- 
tion on the distribution of the data and is a true test 
of modality. The dip probability of the LLS sample be- 
ing bimodal is 88% (91% if limits are removed from the 
sample); hence the distribution is very unlikely to be a 
skewed unimodal distribution. We now apply the GIVLVI 
algorithm to the LLS sample assuming that limits are 
actual values. For the description of the GIVIIVI method, 
its advantages and limitations, we r efer t he reader to 
the appendix in jlVIuratov fc Gnedin| ( |2010[ ). The GMM 
algorithm rejects the unimodal distribution at a con- 
fidence level of better than 99.7%. A heteroscedastic 
split gives /ii = -1.59 ± 0.07, cri = 0.20 ± 0.05 and 
fi2 = -0.38 ± 0.12, (72 = 0.35 ± 0.09 {fii,ai being the 
mean and dispersion of the metallicity, respectively, in 
each branch of the bimodal distribution) . A homoscedas- 
tic spht gives /ii — —1.55 ± 0.08 and ^2 = —0.33 ± 0.10, 
fi = (T2 = 0.28 ± 0.05. The separation of the means 
relative to their widths, D = — ^2|/[(cti -I- a|)/2]°-^, 
is D = 4.5 ± 0.8, which is significantly larger than 2, 
the mi nimum value for a c lean separation between the 
modes (Ashman et al. 1994 1. The two peaks of the distri- 
bution are therefore markedly separated. In both cases, 
the sample is evenly split between the metal-poor and 
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Figure 4. Distribution of the H I column density for the entire 
sample (LLS, SLLS, DLAs). 
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Figure 5. The metallicity of the LLS, SLLS, and DLAs as a func- 
tion of redshift. For the LLS, X is an o element. For the SLLS 
and DLAs, most of the mctallicities were estimated using Zn II. 
The open squares show the [Fe Il/H I] values (error bars are typi- 
cally around ±0.1 dex) for the absorbers where Zn II was not de- 
tected and TVpc u was estimated. The dotted and dashed horizontal 
lines show the average metallicity of the metal-rich and metal-poor 
branches, respectively. Note that metal-poor and metal-rich LLS 
are found at any z at 2 < 1. 

metal-rich groups. If limits are removed from the sam- 
ple, the conclusions would be essentially the same, with 
D = 4.5 ± 0.9, except that the metal-poor group would 
consist of about 36% of the entire sample. Hence the 
GMM method and the dip test demonstrate with a high 
level of confidence that the (unbinned) LLS metallicity 
distribution is bimodal. Note that, as above, fii should 
be considered an upper limit. 

In order to compare the LLS metallicity distribution 
with that of the optically thick absorbers (SLLS and 
DLAs), we systematically searched the literature for 
SLLS and DLAs at 2; < 1 for which the metallicity has 
been estimated. The entire sample is summarized in Ta- 
ble [T] where the absorbers are ordered by increasing A^hi 
values. SLLS and DLAs were largely pre-selected based 
on the observations of strong Mg II AA2796, 2802 ab- 
sorption, with the exce ption of the SLLS and DLAs from 
Tripp et al.| ( |2005[ ) and [Battisti et aH ( |2012[ ). The typi- 
cal selection lor the DLA sample used a criterion of rest 
equivalent width of Mg II A279 6 Wx > 6 A (e .g. 
|et al.„2006j). For this cut-off ,Rao et al.| (|2006|) argHeS" 



might not be the case for the SLLS. We will therefore 
treat the metal-selected SLLS sample with that caveat 
in mind. 

In Fig. |4j we show the H I column density distribu- 
tion for the entire sample summarized in Table [T] There 
are 29 LLS (28 H I-selected, 1 Mg Il-selected), 29 SLLS, 
and 26 DLAs. For the LLS, most of the absorbers have 
16.2 < logA^Hi < 17, while 5/29 have higher A^hi- The 
lack of LLS with log Ah i > 17.5 (tll > 1) is primarily 
an observational bias owing to the UV selection of the 
QSOs, and, secondarily due to the difficulty to reliably 
estimate A'hi at 17.5 < logA'^ni 18.7. This observa- 
tional selection effect therefore affects the H I distribu- 
tion, but there is not an obvious bias imparted to the 
metallicity distribution. 

In Fig. 1^ we show the metallicity of the LLS, SLLS, 
and DLA^as a function of redshift over < z < 1 (i.e., 
over the last 8 Gyr of the cosmic time) . There is no sign 
of evolution within this timeframe (see Fig. [5]). Previous 
authors have shown a general increase of the metallicity 
of the DLAs from z > 5 to z 0, but the metallicity ap- 
pears to plateau from z ~ 1 to z ~ as observed in Fig.p 
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that the dLA sample is essentially complete and there- 
fore free of metallicity bias. On the other hand, this 



DLAs is narrower, with a central value near the present- 
day SMC metaUicity (see next). This suggests the DLAs 
and at least part of the absorbers with logA^ni ^ 20 
have different origins. For the SLLS, the metallicity dis- 
tribution is broader than for the DLAs, with most of the 
values in the range —0.6 < [X/HJslls ^ +0.4. How- 
ever, the distribution might be biased toward the more 
metal-rich absorbers owing to the selection methodology 
for these systems. 

In Fig. [6] we compare the distribution of mctallicities 
for the LLS, SLLS, and DLAs. The DLA metalUcity 
distribution is visually consistent with an unimodal dis- 
tribution with ([X/H]dla) = -0.62 ± 0.35 (1 a disper- 
sion) and the median —0.58. A K-S test rejects the null 
hypothesis of no difference between the DLA and LLS 
datasets at the 99.1% level. The K-S test also informs 
us that the DLA distribution is consistent with a nor- 
mal distribution. In Fig. [7j we have combined the LLS 
and SLLS and separately highlight the absorbers initially 
selected based on their H I content. Metal-selected ab- 
sorbers are unmistakably biased against the low metal- 
licity population; the only detection of metal-poor ab- 
sorbers based on their metal content are strong H I ab- 
sorbers (log Ah I = 18.3 and 20.2). If we consider the 
H I-selected LLS+SLLS, the GMM method would yield 
essentially the same results as listed above; in that case 

* The metallicity reported for the DLAs and SLLS is gen erally 
estimated using Zn II and Fe II (excluding the results from |Bat-| 
|tisti et a l. 2012 and Tripp et al. 2005 where a elements and Fe were 
usedj. As Fe can be depleted significantly onto dust grains and Zn 
is not, we always reported the metallicities derived from Zn. Sev- 
eral SLLS and DLAs have, however, only upper limits quoted for 
Zn based on its non-detection. For these, the abundance deter- 
mined from Fe II does not suggest extremely low metalli cities as 
observed in the LLS, i.e., typically [Fe Il/H I] > —1.3 (see |Meiring| 
|et al.|2009[ |. The studies of [S/Zn]' as a function of the metallicity in 
stars ( non-LTE an alysis) and DLAs show no trend with the metal- 
licity l |Nissen et a l. 2007 Rafelski ct al. 2012). There might be a 
scatter ot about ±0.2 dex, although within la error, [S/Zn] ~ 0, 
i.e., we can directly compare the metallicity distributions derived 
from Zn and a elements. 
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Figure 6. Metallicity distribution function of the LLS, SLLS, 
and DLAs at z < 1. The hashed histograms highlight values that 
are upper limits. The SLLS sample was mostly pre-selected based 
on the presence of strong Mg II absorption and therefore might 
be biased toward the higher metallicity. Most of the DLAs were 
also pre-selected based on strong Mg II absorption, but this should 
not bias as much the metallicity distribution function owing to the 
much larger overall column density. The mean of each peak of the 
LLS metallicity distribution is shown by the vertical dotted lines. 
For the LLS with [X/H] < —1, the mean metallicity is strictly an 
upper limit as 6/14 of the [X/H] estimates are upper limits. 
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Figure 7. Distribution of the metallicity for the LLS and SLLS. 
In this figure, the LLS and SLLS have been combined in one single 
population. The hashed histograms highlight values that are upper 
limits for the LLS+SLLS. This shows that the metal-selected ab- 
sorbers generally miss the metal-poor ([X/H] < —1) branch. The 
two Mg Il-selected absorbers in the metal-poor branch have a large 



H I column density (A^jj i = 10 



J, 10 



19.55 , 



46% of the absorbers are metal-poor. The dip proba- 
biUty of the LLS samples being bimodal is even larger, 
99%. On the other hand, if the entire LLS-fSLLS sample 
is considered, the dip test remains inconclusive. 



0.0 
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Figure 8. The evolution of [C/«] as a function of the metallicity 
of the LLS at 2 < 1 is shown with the circle symbols. The solid 
area shows the approximate t rend and spread of [C/g] vs. [g/H] 
for th e Milky Way stars (from'Akcrman ct al.'2004| [Fabbian et al.| 
|2009| , while the hatched area shows the trend tor z > '2, DLA and 
SLLS (from |Pettini et al.|2008|lPenprase et aLpOlQ^ . The width of 
the areas correspond to about the icr spread at a given metallicity. 



3.2. Relative Abundances of C/a 

We noted in p.3.1| that the relative abundances of O, 
Mg, Si did not depart from solar values within the er- 
rors, and hence, since these elements have a similar nu- 
cleosynthesis history, there is no evidence for significant 
dust depletion (< 0.2 dex) in the LLS. Similar to these 
elements, carbon is not expected to be strongly depleted 
into dust grains at these densities, but the C/a ratio is 
sensitive to nucleosynthesis effects with a time-lag be- 
tween the production of a element s and carbon (e.g., 
Cescutti et aL]|2009| |]VIattsson||2010[ ). In Fig. M we show 
|C/a|-|a;/ii | depenaency determmed from stellar (Ak 



;er 

z > 2 SLLS an 



man et al.|2004 Fabbian ct al. 2009) and 
DLAs (|mtmi~eTral. 2008; Fenp rase et al.||2010[ ) obser 



vations! Chemical evolution models can reproduce that 
trend for [a/H] > —1, but have had mo re problem for 
the upturn of [C/a] a t [a/H] < — 1 (e.g., iCescutti et al. 



2009| )Matts"soii||2010[ ) . Recently, iKobayashTet al.| ( |2011| ) 

show that the observed high [C/a] values in very met ai- 
poor environments are consistent with expected nucle- 
osynthesis yields in core-collapse type-Ill supernova of 
primordial 20-50 Mq stars. 

In Fig.[8l we also show the results for the LLS at z < 1 
where we were able to reliably estimate [C/a]. For the 
LLS with -1.1 < [a/H] < -f 0.1, [C/a] follows a similar 
pattern as the stars i n our Milky Way or H II regions in 
nearby galaxies (e.g., Henry et al. [2000 ). However, for 
[a/H] < -1.4, we find -U.2 < [C/a] < for LLS, higher 
than seen in stars or high redshift DLAs. If there is no 
subtle ionization-dust depletion efi^ect, the solar values 
of [C/a] at metallicities [a/H] < — 1.4 dex could be in- 
dicative of a different chemical evolution history for the 
metal-poor LLS as it scatters around a solar value rather 
than —0.4 to —0.2 dex solar (although we note that all 
the various samples overlap within the observed scatter 
and several values of [a/H] for the LLS are upper limits 
at [a/H] < — 1). The metallicities of the z < 1 LLS are 
so low for these redshifts that most of their metals are 
unlikely to have been incorporated recently into a galaxy, 
although the metals must have been produced and sub- 
sequently ejected from some galaxies at a much earlier 
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epoch. The solar C/a at low metallicity suggests that 
the gas has only been enriched by massive and metal- 
poor stars at an early epoch (see above). 

3.3. Physical Properties 

We now review the physical properties of the LLS and, 
in particular, assess whether the physical properties are 
similar (or not) between the metal-poor and metal-rich 
LLS samples. B ased on the Cloudy models (see Ap- 
pendix and Table Al I, typically the cool gas of the LLS 
is predominantly ionized (/hi: — Nhii/Nh > 90%, and 
on average /hit = 0.98 ±0.02). This contrasts markedly 



from the D LAs that are predominantly neutral (Vladilo 
et al. 2001). From the amount of total iVn (estimated 
trom the Cloudy simulations, see Appendix), we can 
make a rough estimate of the mass of the photoionized 
("phot") CGM gas probed by the strong H I absorbers 
with 16.2 < log Nh I S 18.5. Assuming a spherical 

halo of radius R, its mass is M^^^ w ttR^ fcfimnNn, 
where is the covering factor, /i — 1.3 corrects for 
the presence of He, rriH is the hydrogen mass, and 
A^H = A'^Hi + A^Hii(~ A^Hii for the LLS). We can then 
write: 



M, 



phot 
CGM 



2 X IO^Mq 



R 



0.3 VlSOkpc 
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19, 



From the recent COS-Halos survey of absor bers around 
0.2 < LjL* <i galaxies within p < 150 kpc ( |Werk et ah 
20131, we can estimate the covering factor ot LLS with 
T^IVhi > 16.2: ~ 30% for R < 150 kpc (see also 
Q. The mean impact parameter to galaxies for these 
ciDSorbers is (p) = 60 ± 36 kpc. Based on the results pre- 
sented in the Appendix, for the metal-poor sample the 
logiVn values average around 19.0 ± 0.5, while for the 
more metal-rich sample it is 18.6 ±0.7, implying that the 
metal-rich gas CGM on average may be somewhat less 
massive by a factor ~2 if the sizes and covering factors 
are similar. Hence the gas that is contained in the pho- 
toionized phase traced principally by the strong H I ab- 
sorbers is of 10^-10^° M0. We emphasize that this does 
not include the mass included in the gas-phase traced by 
the observed higher ionization stages or by the lower iVu i 
components at higher absolute velocities relative to the 
strongest A^h i (see below) , which is likely to be at least 
as lar ge if no t much larger ( Tripp et al.|201l" Tumlinson 



et al. 



2011al 



l^'rom the Cloudy models, we can also derive rough esti- 
mates of the linear scale {I = iVn/^-H) of these absorbers 
(see Appendix and references in TablefTl). They generally 
range from a few pc to a few kpc, with one that might 
be larger than 10 kpc. They are therefore not hundreds 
of kpc or Mpc-scale structures as some of t he Ly a forest 



Prochaska et 



ab sorbers (i.e., with lower A^hi. see, e.g., ^ 

al] |2004[ [Lehner et al.| |2006'). There is no trend between 
I and the metallicity (see also Appendix). 

The Cloudy model-derived temperatures are typical of 
photoionized gas temperatures, varying from 0.6 x 10^ K 
to 3x 10"* K (see Appendix and references in TablefTl). As 
the Mg II was observed at high resolution, we can aerive 
reliable 6- values and hence sensitively constrain the tem- 
peratures. We find 3.0 < b < 6.3 kms~^ with the average 
(5) = 4± 1 kms~^ (based on 8 measurements), implying 



T < (0.6-3) X 10^ K. The model-derived temperatures 
are therefore consistent with the observations. 
As we can see from Figs. 



A2 



to 



A18 hot gas (traced 
in particular by O VI, but also by S V, S VI, etc.)|^ is 
often detected at velocities overlapping with low-ion ab- 
sorption in these LLS. We emphasize that although the 
profiles of the low and high ions overlap in velocities, this 
is not sufficient to conclude that the hot and cool ionized 
plasmas have the same origin or meta llicity (see, e.g 



Tripp et al. 2011 Meiring et al.| 2013). To determine 



the metallicity, JSui is required but is very difficult to 
determine for the hot gas. An interesting (but rare) ex- 
ample where the properties of the broad Ly a absorption 
associated with the O VI a bsorption could be es timated 
in a LLS was described in Savage et al.| ( |2011 ). That 
is the absorber at z = 0.2261 toward HE0teT?520 (see 
Table [l]). They show that the hot (- 10^ K) gas traces 
a hot metal-enriched plasma, while the cooler photoion- 
ized gas traced by the LLS is much more metal poor, 
implying different origins (e.g., cool inflow of metal-poor 
gas through a hot metal- enriched halo gas). 

We also note that the presence of hot gas is also not 
ubiquitous along our sightlines. In the bottom panel of 
Fig. [2j we show the total column density of O VI as a 
function of the metallicity of the photoionized gas. Non- 
detections of O VI (down to a level of < lO^'^'^ cm~^) are 
found at any metallicity (25% LLS have no O VI detec- 
tion; 11% LLS have no O VI with hmits < 10"-^ cm'^). 
The scatter seen in this figure may be expected as the 
amount of O VI does not only depend on the abundance 
of O, but of course also on the ionization conditions and 
total hot gas column. For the LLS with A^ovi < lO"'^ 
cm~^. Fig. [2] shows a scatter plot. However, while the 
sample is small, this figure also indicates that LLS with 



O VI 



> 10 



14.5 



cm 



Tiimlinson et al. 

there 



are only found at high metallic 
ity. Such high Nqyi values have been associated with 
galaxies with high star for mation rates and may trace 
the cor onal halo of gala xies ([Heckman et al. "20021 |Chen 

fc Proch aska 2000; Leh ner et a. 

2011a~ Fox||2011( ). Other than this dependence 
does not seem to be a clear dependence between the high 
and low ions for the metal-poor and metal-rich LLS sam- 
ples. A future detailed analysis of the absorption profiles 
of the high and low ions as well as a comparison with the 
output of coUisionally ionized gas models will be neces- 
sary to undertake to assess the relationship between the 
hot and cool gas (or lack of connection) . 

3.4. The Connection to Galaxies 

For some of the LLS in the literature, galaxy searches 
have identified the probable absorber-galaxy association. 
In Table [3j we tabulate in order of increasing LLS metal- 
licity the properties of the LLS-galaxy association for 
9 systems. We only list one galaxy, but there is evi- 
dence toward several sightlines of galaxy interactions or 

^ We have not modeled the highly ionized gas, and therefore 
it might seem unjustified to associate the hig hly ionized gas with 
hot gas. However, based on other works (e.g., 'Lehner ct al. 20091 
Savage et al. 2010, 2011 Tripp ct al. 2011 Tumhnson ct al. 2Ulira 
and on the systematic failure ot our (Jloudy photoionization models 
to reproduce the observed column densities of the highly ionized 
traced by, e.g., O VI, S IV, S V, S VI, we conclude that the highly 
ionized gas found at similar velocities as the LLS must trace coUi- 
sionally ionized gas, likely at T > 10^ K. 
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at least other galaxies ( e.g., 'Jenkin s et al" 2005 Lehner 
et al.|2009 Tumhnson et al..,2011bfe or a possible group 



ot gal axies (e.g., Chen & Frochaska 2000 ; Cooksey et al. 
20081 near the sightlines. In I'ig. [9j we show iVni as 



a function of the impact parameter with respect to the 
associated galaxy for all classes of absorbers, confirm- 
ing that LLS are typically within 100 kpc of a galaxy 
and that stronger A^hi absorbers probe the closer envi- 
ronments of galaxies. In that figure, we did not include 
the z = 0.2823 LLS toward PG1216-h069. This is the 
only absorber for which no galaxy was found within 100 
kpc. In fact the closest galaxy is a 1.85L* at p 3 Mp c 



based on the galaxy survey by Froch aska et al.| ( |2011a 
This galaxy survey is magnitude limited (-K = 19!5), cor 
responding to about > 0.3L* at z ^ 0.3. Hence it is 
quite possible that a less luminous galaxy is closer to 
FG1216-I-069. This highlights a major drawback of the 
current sample of galaxy surveys listed in Tables [l] and 
[3j they are quite heterogeneous in terms of both depth 
and completeness. 

Understanding the galaxy environment is beyond the 
scope of this paper and the issue of galaxy redshift depth 
and completeness will be studied in detail in the fu- 
ture. Nevertheless, despite the current heterogeneous 
and small sample of galaxy surveys, combining our re- 
sults with previous ones, we can conclude that the LLS 
typically probe the CGM of galaxies. The p or Sv param- 
eters do not appear to be a function of the metallicity, 
except possibly that there is more spread in Sv in the 
metal-poor branch (assuming the absorbers are not asso- 
ciated with fainter galaxies). Both sub-i* and super-L* 
galaxies are associated with the LLS in both branches of 
the metallicity distribution, but, again, a better under- 
standing of the completeness will need to be addressed 
to be able to draw stronger conclusions. 

Table 3 

LLS and Possible Associated Galaxies" 



Name 




[X/H] 


P 


Sv 


L 








(kpch-i) 


(kms-i) 


{L*) 


J0943-I-0531 


0.3542 


< -1.30 


95 


+285 


0.4 


PG1630+377 


0.2740 


-1.71 ± 0.06 


37 


-26 


0.3 


TON153 


0.6610 


-1.69 ±0.37 


104 





1.3 


PG1216-(-069'' 


0.2823 


< -1.65 


3238 


-152 


1.9 


PKS1302-102 


0.0985 


< -1.60 


65 


-354 


0.2 


PKS0312-77 


0.2026 


-0.60 ± 0.15 


38 


+16 


0.7 


PG1116-I-215 


0.1385 


-0.50 : 


127 


+50 


2.8 


J1009-I-0713 


0.3588 


-0.40 ± 0.20 


43 


+25 


0.3 


PHL1811 


0.0810 


-0.19 ±0.08 


34 


+36 


0.5 


PKS0405-123 


0.1672 


+0.10 : 


100 


-15 


3.4 


PG1206-I-459 


0.9270 


+0.30 : 


68 





1.8 



Note. — "^We emphasize that the galaxy redshift surveys 
used to find the associated galaxies are heterogeneous in terms 
of both depth and completeness (references for the galaxy sur- 
veys can be found in Table [ill. The reported associated galax- 
ies are generally the closest ro the absorbers in physical space 
(smallest impact parameter p and in redshift spaces (smallest 
Sv = c(2abs ~ %al)/(l + ^abs))- To estimate 5v, the strongest 
H I component is considered. In some cases, several absorp- 
tion components can spread over generally over about 100—200 
kms~^ (and in one case over 1000 kms~^ for the absorber toward 
PG1206+459, see Tripp et al. 2011). We express the abundances 
with the usual logarithmic notation. The references for all the 
sightlines can be found in Table [l] When no information was 
provided on the galaxy luminosity m the published literature, we 
used the information provided to give a rough estimate of L using 
a K-correction. ''Based on a survey complete to about > 0.3L* at 
z ~ 0.3 (see text for more details). 




60 

p (kpc h-') 



Figure 9. The H I column density versus the impact parameter. 
Only the "closest" galaxy (with a spectroscopic redshift) t o a given 
absorber was considered (see Tablelllfor the references and |Meiring| 
|et al.|2011| ). 



4. THE CIRCUMGALACTIC MEDIUM AT z < 1 

For the first time, we are able to probe sensitively 
the metallicity distribution function of the cool (^ lO"' 
K), photoionized (iVni/A^Hii < 0.1) CGM about z < 1 
galaxies (for 16.2 < logA^Hi ^ 18.5). In our analysis we 
only consider the gas associated with the strongest H I 
components, and therefore we have not yet character- 
ized in a systematic manner the connection between the 
strongest and weaker H I components. The weaker H I 
components generally correspond to more highly ionized 
gas seen at higher absolute velocities in the rest frame 
of the LLS (e.g., > 75% LLS have O VI detected at sim- 
ilar or within < 100 kms~^). Analyses of the velocity 
components associated with the LLS show that a large 
fraction of baryonic mass can be present in the lower H I 
column density components where the gas is more ion- 
ized and often highly io nized (e. g., 'Jenkins et al."2005 



Cooksey et aL[2 008; Lehn er etal..2 009; Thom ct al. 201 fl 
Tripp et al."2011; Tumhnso n et al.||2011bj |Kacprzak et] 



ial. „2012b( |Werk et al. 2013 



We find that very low metallicity LLS are not rare at 
z < 1. There is a statistically significant bimodal distri- 
bution of the metallicity of the LLS with a peak around 
[X/H] < —1.57 dex (hereafter the "metal-poor branch") 
and another peak at —0.33 dex (hereafter the "metal-rich 
branch"). Half of the LLS population is found in each 
metallicity branch. However, we point out that metal- 
poor gas found in the same CGM as metal-rich gas along 
the sightline can be effectively hidden if they have simi- 
lar line-of-sight velocities; this effect is illustrated in syn- 
thetic spectra depicting outfiow and inflow in cosmologi- 
cal simulations (Shen et al.,2013) . In contrast to the LLS, 
the metallicities ot the ULAs are normally distributed, 
peaking around [X/H] = —0.6. Not only are the metal- 
licity distribution functions of the LLS and DLAs very 
different, but so is the full range of metallicity values: for 
the DLAs, [X/H] varies from —1.4 to about solar (a fac- 
tor ^ 25 between the lowest and highest values), while 
for the LLS, it ranges from [X/H] < —2 to about -1-0.3 
dex solar (a factor > 200 between the lowest and high- 
est values). The large spread in the LLS metallicities is 
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observed over the entire redshift range < z < 1 (see 
Fig. [7| . The fact that the metallicity distributions of the 
DLAs and LLS are so different implies that they must 
probe different environments or phenomena. In view of 
their larger neutral fraction and smaller impact parame- 
ter to galaxies, DLAs must probe the environments very 
near galaxies, if not the interstellar medium of galaxies 
themselves. 

The bimodal metallicity distribution of the LLS im- 
plies that the gas of the metal-poor and metal-rich LLS 
samples may not appreciably mix in the CGM or, when 
they do, one of the branches must dominate. Indeed if 
some metal-poor gas of mass Mp and metallicity Zp is 
mixed with a metal rich gas of mass Mr and metallicity 
Zr, the metallicity of the mixed gas can be expressed as 

Zrn = (ZpMp + ZrMr)/{Mp + Mr). If Mr » Mp (or 

Mr Mp), then Zm ^ Zr (or Z,„ ^ Zp). If, however, 
Mr ~ Mp, then the metallicity of the mixed gas would 
be somewhat intermediate between the metal-poor and 
metal-rich branches, which is not observed (at least at a 
significant level in the present LLS sample). As we dis- 
cuss below, the metal-poor and metal-rich groups may 
be principally signatures of inflows and outflows, respec- 
tively. In that case, the gas traced by the metal-poor 
and metal-rich branches may not be expected to interact 
much in the CGM. Indeed, in cosmological simulations, 
metal-poor cold flows not only feed the galaxies but also 
supply their angular momentum, and hence they natu- 
rally rotate with the galaxy disks (e.g., D ckcl & Birn- 
boiml2006[ [Stewart et al.|2011bl [Rubin et''al.|2012j |She^ 
■ ' ■■ fe ■■ 



et al. 2013 ) , while on the other hand, stellar feedback in 



galaxies tends to be bi- conical, near perpendicular to the 
galaxy major axis (e.g., Strickland et al.|2004 Bouche et 
aL][2012{ [Brook et al.||20fT| : 

It is ditticult to assess a priori which population of 
galaxies our LLS observations probe. The rate of in- 
cidence of the LLS along the sightline is proportional to 
the product of the comoving number density of galaxies 
giving rise to LLS, tills j and the average physical cross- 
sectio n of the galaxies, ctlls, i-e., dNiA,s/dz oc riLLSO-LLS 



(e.g., Ribaudo et al. 



2011a I. Hence another conse- 



quence of the bimodal metallicity distribution is that 
the product jt-llsO'lls must be similar in both metallic- 
ity branches as about the same number of absorbers are 
observed in each category. While it is well-known that 
the number density of galaxies is related to the galaxy 
luminosity (through the luminosity function), the cross 
section for LLS absorption may also depend on the lumi- 
nosity (e.g. , through a Holmberg s caling, see, e.g., Chen 
et al.[[2010| [Ribaudo et al.[[2011a|). While the number 



density mcreases with decreasmg luminosity, the cross- 
section may decrease. It is a priori unclear which of these 
quantities matters the most for understanding the bi- 
modal metallicity distribution, i.e., whether there could 
be a difference in galaxy type probed by the two branches 
of the LLS metallicity distribution function. 

However, for 10 out of 29 LLS, we know that at least 
one galaxy (sometimes more than one) is found within 
< 100 kpc with a redshift near the LLS (see Tables |T] 
and Is] and Fig. l9|. As we argue in ^3.4 based on previ- 
ous studies, it islikcly that all the LEs probe the CGM 
of galaxies. The metal-rich LLS branch is therefore ev- 
idence that metal enriched gas can be found at several 
tens of kpc from galaxies. The [C/a] ratio follows a pat- 



tern of chemical enrichment that is similar to that ob- 
served in our Milky Way or other nearby galaxies (see 
j |3.2[ ). The metal- rich LLS gas may therefore trace out- 
flows, recycled material, or tidally stripped gas. Detailed 
examples of absorb ers probing outflows can b e found in 
Tripp et al"] \20ll\ or|Prochaska et al.' (2004), while ex- 



ampies o t absorbers tracing g alaxy inte ractions are dis- 
cussed in Jenkins et al. ( 2005 ) and [Lehner et al.[ ([2009 ) . 

The metal-poor LLS population is evidence that near- 
pristine gas can also commonly be found within a few 
tens of kpc from galaxies. The A^hi values of the LLS 
are also too small to trace directly galaxies, and the bulk 
of the metal-poor sample is at [X/H] < —1.4 (only one 
LLS has [X/H] ~ —1-1), i.e., at metallicities too low to 
trace even very low-rnetallicity galaxies at z < 1 (e.g . , 
Kobulnicky fc Kewley|[2004[ [Perez- Montero et a.l.[|2013[) , 
with only exceptions being too uncommon to find with 
regularity in absorption line studies (see also 
|aL^2011bj). 



Ribaudo et 

The finding of low metallicity LLb with near 



solar C/a (see ^ 3.2 1 is also consistent with the metal-poor 



CGM gas being unlikely to have been recently incorpo- 
rated into galaxies. Thus they may instead trace cold 
streams of ionized, metal-poor gas falling onto galaxies 
with properties consistent with those predicted in cos- 
mological simulations (e.g., Fumagalli et al.[[2011b Shen[ 
let al.,,2013. see below ). 



Werk et al.[ (2013) have reported some properties (but 
not the metallicity) of the absorbers within 150 kpc of 
42 galaxies from the COS-Halos project, probing a range 
of luminosity 0.2 < L/L* < 3 (about 60% of galaxies be- 
ing < 0.7 L*). The COS-Halos sample was not selected 
based on the A'hi content, but on the galaxy proper- 
ties and impact parameter. The A^hi distribution of this 
sample follows more clos ely the differential H I column 



density distribution (e.g., Lehner et al.|2007 Ribaudo et 



al. 2011a) , and about 70% (including tne lower limits) ot 
their absorbers have 14.5 < log A^hi < 16- About 40% of 
those have no detection of Mg II, but in this A'h i inter- 
val, they are not sensitive to metallicities [X/H] < — 1. 
About 30% of the absorbers have logA^ni > 16, but only 
lower limits on A'hi have been determined . The only 
exception is one absorber, which is in our sarnp le, the 
metal-poor LLS toward J0943-I-0531 (see Table [l} iThom 
et al. 2011). It is therefore too early to determine if there 
is any difference between the COS-Halos and our samples 
regarding the metallicity distribution of the LLS. Fur- 
thermore since that there are only about 10 COS-Halos 
absorbers that are in the LLS A^hi interval, a statisti- 
cally robust comparison between the two samples may 
actually be impossible to undertake. 

Other observational techniques have also been used 
to empirically characterize the flows of matter between 
the galaxies and the CGM. One, in particular, has 
been extremely useful for characterizing the outflows of 
galaxies, whereby "down the barrel" spectra of galaxies 
are searched for blueshifted and redshifted absorption 
troughs relative to the galaxy systemic velocity as signa- 
tures of outflow and inflow, respectively (e.g. , Hcckman] 
let al.''2000; ■Martin"2005'; 'Tremonti et al."2007; WcineTj 
iet al. 2009; S tcidel et al 2010 ; Rubin et al. 2011). This 
methdcT reveliis that outflows are quite ubiquitous in star- 
forming galaxies. However, with this technique and over 
the redshift interval 0.4 < z < 1.4 that overlaps w ith our 



study, Rubin et al. (2012) and Martin et al. (2012) found 
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a very small fraction (only a few percent) of galaxies 
showing significant redshifted (infalling) Fe II or Mg II 
absorption relative to the host galaxyp"] As these au- 
thors argue, the gas that is seen is metal enriched, and 
hence more likely tracing the recycling of gas via, e.g., 
a galactic fountain. Identifying infalling metal-poor gas 
has proven difficult using metal-lines (especially Mg II 
and Fe II). Our results demonstrate there is a significant 
amount of this metal-poor gas, but that it would produce 
Fe II and Mg II absorption too weak to be identified in 
complicated galaxy spectra where outflows are very often 
present. While the QSO absorption line technique does 
not provide the direction of the gas motion relative to 
the associated galaxy, and requires additional observa- 
tions to characterize the galaxy properties, selecting ab- 
sorbers drawn from an H I-selected population is unique 
in providing insights into the interaction of galaxies with 
their surroundings and may be the only way to directly 
study th e very low metallicit y CGM in absorption (see, 
however, Ranch et al.|[2011| for Lya emission imaging 
t echniques^ 



the equations fSpringel'12010 Vogelsberger et al. 2012) 
These simulations show m particular some critical dif 
ferences in the covering fac tors of LLS and DLA and 



The present work and other recent r esults (e.g., Stocke 
eraL||20T0l [20T3l |Thom et al.|[20TTl I'l^ipp et al ||l^011 

lLI 



Tumlinson et alJ2 011a[ |Ka cprzak et al.|2012b[ Church! 
et al. 2012; 'Werk et al. 20l3p signal a snift in our ability 
to characterize the properties of the CGM and have al- 
ready provided new stringent empirical results to test 
cosmological hydrodynamical simulations. Regarding 
the metallicity of the CGM, there are now several theo- 
retical works that determine the metallicity of the inflow- 
ing and outflowing gas at z ~ 2-3, although there is not 
yet a theoretical prediction of the metallicity distribution 
function of the CGM in the LLS regime at z < 1 for com- 
pa rison with the present re sults. The recent simulations 
by Fumagalli et al. (2011b I include mild stellar feedback 
and show that cold streams are traced mostly by LLS 
within 1 or 2 virial radii of galaxies, where the gas has 
only been weakly enriched ([X/H] ~ — 1.8±0.5). The gas 
is also pred ominantly ionized a s it is observed. The sim- 
ulations by Shen et al. (2013) include stronger galactic 
outflows, but still show that cold flows are metal-poor 
with a median value —1.2 dex. These simulations and 



those by van de Voort & Schaye (2012) show that major 
outflows do not suppress signiti cantiy cold flow accretion, 
van de Voort & Schaye (2012) similarly show that cold 
mode accretion is generally rnetal-poor ([X/H] < —1.5) 
for any halo mass. Our results for the metal-rich LLS at 
z < 1 show metallicities higher than predicted in their 
models, but their simulations are at z = 2. The metal- 
licities of the cold flow gas from these theoretical results 
compare quite well with our empirical results, and it is 
unlikely that the metallicity of the cold fllaments is much 
dependent on the redshift. However, it would be im- 
portant to confirm this in the simulations as well as to 
understand the type of galaxies and their properties as- 
sociated with these absorbers within these simulations. 
While this comparison between the simulations and the 
present observational results is encouraging, it will also 
be important to determine if these recent theoretical re- 
sults are confirmed with the AREPO simulations where 
different hydrodynamic schemes are employed to solve 



ho w galaxies get their gas ( Keres et al. 2012 Bird et 



al. 2013 Nelson et al. 2013)! INotwithstanding, and in- 



dejaendently of the theoretical work, our observational 
results show that large amounts of both metal-poor and 
metal-enriched cool, photoionized gas traced by the LLS 
are found around galaxies at z < 1. 

We finally point out that the often adopted 10% solar 
metallicity as a representative value for the IGM in the 
low redshift universe has no support from our observa- 
tions. In fact, the O VI absorbers often used to charac- 
terize the warm-hot ionized medium (WHIM) are quite 
likely to trace the CGM rather th an the IGM based on 



their overall high metallicity (e.g., [Sembach et al. 2004 



Although [Rubin et al.| also noted that the presence of infall 
onto as many as 4U7o ot the galaxies could not be ruled out. 



Lehner et al. 2006 1, consistent with the metal-rich branch 
of the LLS. 'rhe rnajority of the O VI absorbers observed 
at z < 0.5 are also associated with sub-dwarf galaxies, 
not th e diffuse WHIM predict ed by co smological simula- 
tions ([Prochaska et al. 201 lb[ see also Wakker & Savage 
2009 Stocke et al. 2006). Another consequence ot this 
work IS therefore that the IGM far from the virial radii 
of galaxies could have a characteristic metallicity well 
below 2% solar at z < 1. The very low metallicity gas 
has been essentially missed in studies of the Lya forest 
because absorbers with logiV(H I) < 14 have a column 
density too small to detect metal lines at significant levels 
to constrain the metallicity of the gas below 10% solar. 

5. SUMMARY 

We have assembled the first sizable sample of LLS 
at z < 1 selected solely on their H I content (16.2 < 
logA^Hi ^ 18.5), i.e., where the main criterion for as- 
sembling our sample of CGM absorbers is independent 
of metallicity. The iVn i of the LLS is too low to directly 
probe the galaxies, but is large enough to trace their 
CGM. The metallicity is determined by estimating the 
column densities of metal ions and H I in the strongest 
H I component (not over the entire velocity profile where 
metal- line absorption may be observed), and correcting 
for ionization effects using Cloudy photoionization mod- 
els. There are 28 LLS in our sample, 16 of them being 
newly analyzed in this work. Our main results are as 
follows: 

1. We empirically establish the metallicity distribu- 
tion function of the LLS (and hence the cool CGM) 
at z < 1: the distribution is bimodal with a peak 
around [X/H] < -1.57 ± 0.24 (the metal-poor 
branch) and another peak at [X/H] = —0.33 ±0.33 
(the metal-rich branch) with about equal number 
in each branch. The metallicity in the metal-poor 
branch is so low at z < 1 that the majority of the 
gas in in these LLS is unlikely to have been recently 
incorporated into galaxies. 

2. In contrast, the DLA metallicities are normally dis- 
tributed, peaking around [X/H] ~ —0.6. The range 
of metallicity between the DLA and LLS samples 
is also quite different, varying from [X/H] — —1.4 
to about (a factor ^ 25 between the lowest and 
highest values) for the DLAs, while for the LLS, it 
varies from < —2 dex to about -1-0.3 dex (a factor 
> 200 between the lowest and highest values). 



The Bimodal Metallicity Distribution of the CGM 



17 



3. The gas associated with the strongest H I com- 
ponent of the LLS is predominantly photoionized 
(/hii > 90%), with a temperature around 10'' K. 
The mass of the photoionized component of a typ- 
ical galaxy's CGM selected by these LLS is ~ 10®- 
10^" M0. There is also often (but not always) evi- 
dence of more highly ionized gas as revealed by the 
presence of, e.g., O VI, S VI, and other highly ion- 
ized species as well as other higher absolute velocity 
components with lower A'hi, which may contain as 
much mass. 

4. We find that C/a varies with metallicity for the 
z < 1 LLS. For [a/H] > -1, C/a follows the gen- 
eral pattern seen in stars and H II regions with 
similar metallicities, suggesting a similar chemical 
enrichment. For [a/H] < —1, G/a is, however, near 
solar, different from the trend seen in low metallic- 
ity stars and high redshift DLAs. The solar C/a 
at low metallicity suggests that the gas has only 
been enriched by massive and metal-poor stars at 
an early epoch. 

5. The bimodal metallicity distribution of the LLS im- 
plies that the metal-poor and metal-rich samples 
have similar tills clls and that the metal mixing 
between the two metallicity branches is either inef- 
ficient or dominated by one branch. 

6. Combining these results, we conclude that the 
metal-rich gas likely traces winds and recycled gas 
from outflows and galaxy interactions. The metal- 
poor LLS are very likely tracing cold accretion onto 
galaxies, with properties in very good agreement 
with those seen in cold flow accretion models and 
simulations. Independent from the simulations, our 
empirical results show there is not only a large mass 
of metal- rich gas around galaxies at z < 1, but also 
a significant mass of metal-poor gas that may be- 
come available for star formation. 

7. Finally, based on these results, the metallicity of 
the diffuse IGM at z < 1 may be well below 2% 
solar. 
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APPENDIX 

In this Appendix, we provide more details for the new sample of LLS, in particular about their velocity profiles 
and column densities, and the Cloudy photoionization models that were used to determine the metallicity of these 
absorbers. The redshift of each absorber is defined based on the strongest H I component. In the figures that follow, 
we show the normalized profiles of H I (generally a weak and a strong transition) and most of the detectable metal 
lines observed by COS as a function of the rest-frame velocity. We also systematically show the Mg II transitions when 
optical observations were acquired. 

It is useful to first describe the Cloudy models for two illustrative systems, the z — 0.5185 and z = 0.7289 absorbers 



toward PG1522-I-101. The results of Cloudy modeling of these figures are shown in Fig. Al These plots show the 
predictions of Cloudy models for an adopted metallicity as a functio n of ionization parameter U . The Cloudy results 
for these two and other absorbers are summarized in Tables [T] and |A1| for the new sample of LLS described here. 
We also list in these tables the results from the existing sample of LLS at z < 1 prior to this work (the references 
are provided in Table [T]). For each absorber, the column densities determined from the observations and the Cloudy 
simulations are matched by varying the ionization parameter ([/) and metallicity. The errors that we provide on U 
and the metallicity are not formally la error, but are ranges of values allowed by the models and errors on the column 
densities. From the Cloudy solution, several other useful physical parameters can be determined (including, the total 
H column density (iVn), the temperature of the gas - T, the density of the gas - nn, the pressure of the gas - P/k 
and the linear scale of the cloud - I = A^h/"^h)- The values corresponding to the best solution are listed in Table Al 
A description of these two absorbers follows. We comment below on all of the absorbers in our sample, but only snow 
the models graphically for these first two, which were chosen because they are characteristics of the high and low 



metallicity branches o f the metallic i ty dis tributions. Other models can also be found in, e.g., Lehner et al. (2009) 



Ribaudo et al | ( |2011b| ), |Tripp et al (20111. 

- PG152^01, z = 0.5185: I-'ig. |A2| shows that the Mg II absorption has more than one component, but the compo- 
nent at kms""'^ dominates the absorption. This is also the case for H I where most of the absorption in the weakest 
H I transitions is observed at kms^^. For O II, O III, C II, and C III, there is no evidence of multiple components at 
the COS resolution. The strongest H I component aligns well with the strong Mg II component as well as the metal 
ions observed with COS accounting for a COS velocity calibration error of ~5-15 kms-\ A single component COG 
fit to the H I lines and the AOD measurements of the weakest H I lines are consistent within la. We adopt a mean 
from the results of these two methods (in this case, the continuum break is too weak to estimate iVni reliably). 
Based on the com parison of the velocity profiles of the different ions and H I, a single-phase photoionization model 



is justified. Fig. Al shows that a Cloudy photoionization model with \ogU = —3.6 ±0.1 and [X/H] = —0.40 ±0.05 can 
reproduce well the O II/Mg II, C II/O II, O II/O III ionic column density ratios. The well determined C III column 
density suggests that [C/a] ~ —0.15, consistent with Type II supernovae enrichment. As there is no evidence for 
highly ionized gas, this absorber is single phase. 



PG1522+101, z = 0.7289: Fig. A3 shows that there is no strong evidence for multiple components in the H I, C III, 
or O III absorption profiles. There is no detection of any singly ionized species (including the strong Mg II A2796 line) 
and S III. C III A977 is partially contaminated at higher velocities, but that does not affect much the component at 
kms"^ Besides C III, O III A832 and O VI A1031 are the only reliable metal line detections. All the useable H I lines 
(i.e, lines that are uncontaminated and for which the continua can be reliably determined) remain quite strong, and 
therefore the AOD method only provides a lower limit on A^hi- We determine iVfji from the continuum break and 
the COG. Both methods give consistent A^hi values within la (and the AOD limit is consistent with these estimates), 
and we therefore adopt a mean from these two methods. 
Based on the comparison of the velocity profiles of the different ions, we used a si ngle phase Cloudy photoionization 



model to predict the column densities or limits of all the observed species (see Fig. Al ). A model with [X/H] < —2.0 
and \ogU > —3.20 can reproduce the observations. Higher metallicity would overproduce Mg II and S III. Note that 
even though we detect metals (C III and O III), the lower limit for U leaves the total H poorly constrained, giving 
only an upper limit on the metallicity. In this case, C/a is consistent with a solar despite the gas being extremely 
sub-solar. The model cannot produce enough O VI (log A^ovi = 13.92 ±0.07, but we emphasize that this is only based 
on the strong transition of the doublet since O VI A 1037 is contaminated), and hence this absorber is multiphase. As 
we discuss in the main text, for the LLS, photoionization alone fails to reproduce the columns of the high ions when 
they are detected, implying that other processes must be at play (e.g., coUisional ionization). 

For the remaining 14 LLS, similar Cloudy simulations and figures were produced, but we only discuss the results 
below. The order of the absorbers follow that of Tabled 



- PHL1377, z = 0.7930: Fig. A4 shows that Mg II A2796 is barely detected at km s ^ (the weaker transition at 
2803 A is not detected at 3a). Most of the H I col umn density is in the kms ^ component, but there is evidence 
of at least another component around +40 kms~^ in the absorption profiles of C III, O III, O IV, and O VI (Si IV is, 
however, mostly present at kms~^). The COG, AOD, and continuum break methods all give consistent results for 
Nhi and our adopted A'^hi value is a mean of these results. O III, O IV, C III have all large central optical depth, tq, 
suggesting that the saturation might be important, and hence we treat the column densities for these ions as lower 
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limits. C IV A1550 appears slightly contaminated at u > 70 kms , and we therefore used only C IV A1548 to estimate 
Nc IV (note that the C IV data - as weU as Si IV and Al II- were obtained with STIS E230M). O VI A1037 is not 
covered at this redshift, and therefore only O VI A1031 is used to estimate iVovi- 

As the strong doubly ionized species (C III, O III) show at least two strong components, we focused mostly on the 
singly ionized species and weak S III transitions for the ionization model. The total column densities of C II, O II, Mg II, 
and S III can be satisfied simultaneously within la error in the Cloudy photoionization model [X/H] — —1.45 ± 0.05 
and \ogU — — 2.90±0.10. Using, C II, O II, and Mg II, we derive [C/a] ~ —0.05 despite the low metallicity of the gas. 
This model produces enough C III to match the observed lower limit, about 50-60% of the observed O III and Si IV, 
and much too small amounts of O IV, O VI, C IV, and S V relative to the columns determined from the observations. 
This absorber is therefore multiphase. 



PG1338+416, z = 0.3488: Fig. A5 shows that the Mg II absorption profiles reveal two main components at and 



+35 kms ^, but the column density is dominated by the kms ^ component. The H I has at least two components at 
and about —80 kms^^, but 90% of the column density is in the kms~^ component. The one component COG fit 
to the H I lines gives log A^hi = 16.40 ± 0.05, while the AOD results on the weakest measurable H I line (A916) yields 
log A^Hi — 16.20 ± 0.07. The difference could be due to the presence of the +35 kms^ component in the strong H I 
transitions that the COS data does not really resolve. We adopted a mean of the two results, where in this case the 
adopted error reflects the scatter between these two values (no reliable estimate could be made from the continuum 
break owing to the proximity of the Galactic Lya absorption component and the unreliable flux level near the edge 
of the detector). The C III A977 and Si III A1206 profiles have strong absorption and may be saturated (although our 
Cloudy simulation implies the saturation should be mild, < 0.2 dex). The O VI AA1031, 1037 profiles are strong as 
well, but the lines are unsaturated based on the excellent Na agreement between the O VI doublet lines. A similar 
agreement is found for the Mg II doublet. 

The total column densities of Mg II, Si II, and Si III can be satisfied simultaneously in the Cloudy photoionization 
model with [X/H] = -0.75 ± 0.15 and logJ7 = -3.2 ± 0.10. If C II is not contaminated, then [C/Si] ~ +0.15, which 
is higher than typically observed at this metallicity. The photoionization model cannot reproduce Nqyi by orders of 
magnitude, implying that the absorber is multiphase. 

- PG1338+416, z — 0.6863: Fig. A6 shows that the Mg II absorption profiles reveal two main components at and 

45 kms"""^, but the column is dominated by the kms~^ component. This conclusion also applies to the other singly 
ionized species. However, for the doubly and higher ions, the absorption is dominated by higher velocity components 
(e.g., about 45 kms~^ for C III and O III, and 70 kms~^ for O VI, S V, S VI). The H I has at least two components at 
and about 45 kms~ , but most of the column density is in the kms^ component. The single component COG fit 
on the H I lines gives logA'ni = 16.39 + 0.02. The AOD results on the weakest uncontaminated H I lines (AA916, 917, 
919, 920) yields logTVni = 16.45 + 0.04. From the continuum break, we estimate log A^hi = 16.48 + 0.06. We adopted 
a mean for Ahi from these 3 methods. For C III A977, O III A702, and O IV A787, the absorption is strong and likely 
saturated. Based on the difference of 0.07 dex between the weak and strong lines of the O VI and Mg II doublets, we 
applied a 0.07 dex correction on the weakest line of each doublet ( Savage fc Sembach]|1991 ). 

Based on the comparison of the velocity profiles, with the Cloudy simulations we only attempted to reproduce the 
column densities of C II, O II, and Mg II. The total column density of the ions can be satisfied simultaneously in the 
photoionization model with [X/H] = +0.10 ± 0.15 and \ogU = —3.90 ± 0.15. The metallicity cannot be much below 
< —0.1 dex because otherwise not enough Mg II would be produced for any U . If U is higher, then O II would be 
overproduced by +0.5 dex. For that absorber, [C/Si] ~ +0.15, which is typical for super-solar metallicity gas. We 
note that even if U is small, the absorber is still nearly pre-dominantly ionized (more than 90% ionized). The linear 
scale is, however, quite small, of the order of 3 pc. Its small linear scale implies that it is very unlikely a "clump" or 
"cloud" of gas. Instead it might a somewhat dense "shell" within an extended ionized medium (possibly probed by 
the other velocity components). This absorber has more than one component and is multiphase. 

- PG1407+265, z = 0.6828: Although the H I C III, O III, and O IV profiles are strong and suggest two components 

at and about —30 kms^^ (see Fig. A7), the profiles are dominated by the component at kms~^ as evidenced by the 
weaker transitions of H I and S III AA698, 724 lines. The strongest transitions of O II and Mg II are not detected, but 
O VI and S V are also observed at kms^^. The C III absorption is quite strong and could be somewhat saturated. 
The COG, AOD, and continuum break methods all give consistent results for Ahi and our adopted result is a mean 
of these results. All the ions with two transitions (O VI, S III, O III) give consistent column density results. 

For this absorber, there is no evidence for very different kinematics between the different ions. We therefore con- 
sidered all the limits and detections for the Cloudy model. The non-detection of Mg II and C II and the detection of 
S III constrain this absorber to —2.0 < [X/H] < —1.5 with —2.3 < \ogU < —1.7. In the case of the highest value of U 
(lowest metallicity), all the observations except O VI can be matched. For the lowest value of U, the model can only 
reproduce the singly and doubly ionized species. Depending if C III is saturated or not, then —0.2 < [C/alpha] < +0.2. 
O IV, O VI, S V, and S VI, cannot be r eproduced by photoionization, and therefore this absorber is also multiphase. 

- PKS0552-640, z = 0.3451: Fig. XS shows that only the metal lines Si III and C III are detected in this absorber, 

both with only one component at the COS resolution. C III is blended with another feature at w < —30 kms~^, but 
based on Si III we do not expect the absorption to extend much below v < —30 kms~^. The single component COG 
fit of the H I lines gives log A^hi = 16.83 ± 0.04, overlapping within 2a with the results from the continuum break. 



The Bimodal Metallicity Distribution of the CGM 



21 



logiVni = 16.95 ± 0.05. The combination of the hnes being strong (the weakest hnes being on the knee of the COG) 
and the break occurring in the wing of the damped Galactic Ly a hkely explain these differences. We adopted a mean 
for iVni between these 2 methods, where the adopted error reflects the scatter between these two values. 

Based on the comparison of the velocity between the different ions, we used a single-phase Cloudy photoionization 
simulation to model all the column densities of the observed species and limits for the non-detected species. A model 
with [X/H] < —1.60 and logJ7 > —3.70 can reproduce the observations. In this case, [C/a] ~ —0.15. We emphasize 
that this absorber could be an extremely low metallicity absorber, and a better limit on Mg II would be extremely 
beneficial to better constrain the metallicity. As there is no O VI detected or other ions that cannot be produced by 
the photoionization model, it is a single-phase absorber. 

- PKS0637-752, z = 0.4685: The profiles of O II and weak H I transitions suggest that they are dominated by a single 
component (see Fig. A9), while the C III and O III absorption profiles show an additional but comparatively weak 
component at about —50 kms~^. The two detected unblended O II hnes and the O VI doublet lines produce consistent 
column densities for these ions, respectively. The C III and O III lines are relatively strong and may be saturated. The 
Mg II doublet was only observed at i? ~ 5, 500, and there is evidence for weak satura tion comparing the results from 
the weak and strong transitions. We corrected for the saturation (0.05 dex) following Savage & Sembach (1991). The 
single component COG fit on the H I lines gives logA'ni — 16.44 ± 0.06, overlapping within Icr with the result from 
the continuum break, log A^hi = 16.51 ± 0.05. We adopted a mean for A^hi between these 2 methods. 

For this absorber, we can estimate a strong lower limit on the metallicity, [Mg/H] > —0.5, based on Mg II. A model 
with [X/H] = —0.4 ± 0.1 and logC/ = —3.8 ± 0.1 provides a good agreement between the observed and modeled 
column densities of Mg II, O II, C II, and C III. Aq hi is, however, underproduced by a factor 6, possibly because it is 
partially contaminated or because the gas is multiphase as evidenced by the absorption of O IV and O VI (which are 
underproduced in the photoionization model by several orders of magnitude). A higher metallicity (and smaller [/) 
would be consistent with the singl y ioni zed species, but not C III, and we therefore favor the above solution. 

- HE0439-5254, z — 0.6151: Fig. AlO shows that the profiles of this absorber have a complicated velocity structure, 

with the absorption spanning from about —160 to +200 kms~^. However, most of the H I column density is located in 
the component at kms~^ between —75 and -|-70 kms~^ as evidenced by the weakest H I transitions. The weak H I 
transitions also show a weaker component at -1-145 kms~^, which is cleanly separated from the low velocity component. 
The singly ionized species and S III profiles have a similar velocity structure than the weak H I transition. The higher 
ions with strong transitions, including C III, O III, O IV, and O VI, show absorption over velocities comparable to that 
in the Lya absorption profile. The profiles of the high ions with weaker absorption (S IV, S V, S VI) is dominated 
by the component at kms~ . Therefore, except for C III, O III, O IV, and O VI, the column densities of the metals 
and weak H I transitions can be readily compared in the kms~^ component between —75 and -1-70 kms^^. The 
one component COG fit to the H I lines (from A916 to A949; at higher wavelengths all the components are entirely 
blended) and AOD estimates using the weakest H I lines (A < 923) give a similar A'hi. We adopted a mean of the 
results. 

We only attempted in our Cloudy simulations to predict the column densities of C II, O II, and S III where the 
absorption is mostly confined to the strongest H I component at kms~^, where the H I column density is estimated. 
A model with [X/H] = -0.30 ± 0.05, logf/ = -2.70 ± 0.10, and [C/a] = -H0.05 ± 0.10 reproduces very weh the 
observations. In that case, the model predicts too much O III, which is consistent with the lower limit determined by 
the observations. This model, however, only produces a small amount of O IV and S IV and negligible amount of S V. 
It is therefore a multiphase absorber. 

- SBS1122-I-594, Zabs = 0.5574: The velocity profiles for this absorber have a complicated structure (see Fig. All), 

with the absorption spanning from about —200 to -1-120 kms^^. However, most of the H I column density is located in 
the component at kms~^ based on the weak H I transitions. Within this component, there is only a weak absorption 
present in the profiles of the metal lines. The ionization and/or metallicity conditions must be different between the 
various components since a very strong C III and O III absorption (a factor > 6 times larger in column than in the 
kms~ component) is observed at about —125 kms" with a smaller A^hi (log Ahi — 15.85 or factor 2.8 smaller than 
that of the kms~ component). In this case, we can estimate the column densities in the individual components and 
we only consider the component at kms~ , i.e., the absorption between about ±50 kms~ . C II A904 is a tentative 
detection (C II A1036 is not detected, but the limit is consistent with A'cii estimated at 904 A). The one component 
COG fit to the H I lines (from A916 to A937) and AOD estimates using the weakest H I lines (A < 920 A) give the 
same A^hi- We adopted a mean from the two results. 

This absorber is complicated to model because C/0 may not be solar, so we have to rely solely on the O II/O III 
ratio (estimated between —50 and -t-50 kms^^) to estimate the metallicity. While the O II absorption is weak, the 
two available transitions at 833 and 834 A provide a consistent A'oii- The detection of O II provides a strong lower 
limit on the metallicity [X/H] > —1.1. Assuming that all the O III and O II are in the same gas phase, then both 
the metallicity and U are tightly constrained with [X/H] — —1.05 ± 0.05 and logJ7 = —3.2 ± 0.10. For that solution, 
we need to have [C/0] ~ —0.6. O VI may be present at kms^^, and therefore this absorber is multiphase, as the 
observed A^o vi cannot be reproduced by the photoionization model. 
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- J1419+4207, ^abs = 0.2889: While the high resolution Mg II AA2796, 2803 spectra shown in Fig. [Il2| reveal three 

components at —14,0,24 kms~^ (the latter being the weakest), that velocity structure cannot be discerned in the 
lower resolution COS data. We therefore report the total column densities, i.e., Na{v) is integrated over the entire 
velocity profiles for that absorber. C III is quite strong and could be saturated. The one component COG fit to the 
H I lines overlaps within la with the AOD result. We adopted a mean from the two results. 

A Cloudy photoionization model with [X/H] = -0.65 ± 0.10, \ogU = -3.10 ± 0.10, and [G/a] ~ -0.15 reproduces 
simultaneously the column densities of C II, C III, Mg II, and Si III within the estimated la observational errors 
(allowing for about 0.2 dex saturation for C III), and the limit on Si II. A higher metallicity would over-predict the 
amount of Si II. A lower metallicity would under-predict the amount of Mg II. These two constraints set the t/ in a 
narrow range of allowed values. O VI cannot be produced at all by this photoionization model, so it is a multiphase 
absorber. 

- J1419-|-4207,Zabs = 0.4256: No Mg II is detected in this absorber (see Fig. |A13|, and the COS absorption profiles 
do not reveal more than one component. All the metal lines are relatively weak, suggesting no saturation issue. The 
one component COG fit to the H I lines and AOD estimates using the weakest H I lines (A < 920 A) give the same 
A^Hi. We adopted a mean of the two results. 

A Cloudy photoionization model with [X/H] — — 1.40±0.20 and \ogU = — 2.90±0.20 reproduces simultaneously the 
observed column densities of Si III, C III, and agrees with the limit on Mg II. A higher metallicity would over-predict 
the amount of Mg II. A lower metallicity would under-predict the amount of Si III. O VI cannot be produced at all by 
this absorber, so it is a multiphase absorber (however, we note that O VI is only detected in the stronger line, is weak, 
and its peak optical depth is shifted relative to the low ions; all this leading to an uncertain detection of O VI). 

- J1419-I-4207, 2;abs = 0.5346: Fig. A14 shows that a strong single Mg II component is observed in this absorber. 



The only other significant detection of a metal line is C III. The O II AA832, 834 absorption lines are weak, but both 
transitions give consistent results. C II A904 and O III A832 are barely 3a detections. We estimated the same apparent 
column density for the weak and strong Mg II lines, implying that the lines are fully resolved. Unfortunately, for this 
absorber the H I column density remains not very well determined (relative to the other absorbers) because the S /N is 
too low to provide reliable estimates of Na or W\ below 920 A: so all the available lines are near or on the flat part of 
the COG. The S/N is also too low to determine A^hi from the continuum break method. The adopted result is from 
the one component COG fit to the H I lines. 

A Cloudy photoionization model with [X/H] — —0.20 ± 0.20 and log U — —3.9 ± 0.20 reproduces simultaneously the 
column densities of O II, C III, and Mg II. If O III is real, then the modeled column density is about a factor 5 too 
small. This is similar to the absorber aX z — 0.4685 toward PKS0637-752. No O VI is detected (although only down 
to a level of 13.90 dex owing to the l ower S/N spectrum of this QSO), but if O III is real, the gas could be multiphase. 

- J1435-I-3604, Zabs = 0.3730: Fig. A15 shows that a very weak single Mg II component is observed. The only other 
detection of a metal line is C III. The AOD estimate of the weakest measurable H I transition (A916) and the COG fit 
to all the uncontaminated H I lines give the same A^hi- The estimate of the column density break is consistent with 
the cumulative effect of this absorber with the absorber at z = 0.3878. 

Based on the Cloudy photoionization model, the detection of Mg II implies that [X/H] > —2.0 as below that value, 
not enough Mg II is produced for any U. If [X/H] > -1.7, then \ogU < -4.1 and [C/Mg] > -1-1.2, the latter being 
very unlikely. Therefore there is a small range of allowed metallicity that can reproduce the column densities of C III 
and Mg II (and limits on other ions), [X/H] = -1.85 ± 0.10 and logt/ = -2.9 ± 0.20. No O VI (log A^ovi < 13.51) is 
observed in this absorber, so there is no evidence for a multiphase absorber. 

- J1435-f 3604, Zabs = 0.3878: As shown in Fig. A16 no Mg II absorption is detected. The only detected metal lines 
are C III and O III. Both have more than one component, with the optical depth peaking in the second component 
that is shifted by about —40 kms~^. Similarly O VI A1037 is also shifted by the same amount from the strongest H I 
component. The one component COG fit to the H I lines and AOD estimates using the weakest H I lines (A < 930 A) 
yield similar A^hi. We adopted a mean of the two results. 

Based on the Cloudy photoionization model, the non-detections of Mg II and Si III place a strong upper limit on 
the metallicity: for any U, [X/H] < —1.4. This also consistent with the 3a upper limits on O II, although that limit is 
not very constraining. We are unable to better constrain the metallicity because, first, A^oiii ^nd N^m are somewhat 
uncertain over the velocities where A^hi is estimated, and, second, the combination of increasing U and decreasing the 
metallicity can fit all the observables. In order to be able to produce A'oiii and Nqui, logt/ must be > —3, but we 
emphasize that the velocities of these ions are shifted relative to the strongest H I component. This absorber is also 
multiphase. 

- J1619-I-3342, Zabs = 0.2694: Fig. A17 shows that the Mg II absorption is very weak but detected in the stronger 
line of the doublet. There are possibly two components. The absorption of C II A1334 is barely detected at 3a, and 
Si II A1260 and Si IV A1393 are not detected. In contrast, C III and Si III are both well detected and line up well with 
the H I absorption within the COS velocity calibration error. A single component COG fit to the H I lines, the AOD 
measurement of H I AA917, 920, and the continuum break method all give overlapping results for A'hi within la. We 
adopted a mean for A'hi of the results from these 3 methods. 

A Cloudy photoionization model with [X/H] = -1.60 ± 0.10, \ogU = -2.9 ± 0.15, and [C/a] ~ -0.10 reproduces 
simultaneously the column densities of Si III, Mg II, and C III and limits on C II. The metallicity cannot be higher 
as otherwise for any U, too much Mg II would be produced. There might be a 2-3 cr detection of O VI A1037, but 
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unfortunately O VI A1031 is contaminated, therefore it is unclear if this absorber is multiphase. 

- PG1216+069, Zabs = 0.2823: For that absorber, we use both the COS and STIS E140M observations (for a detailed 
description of the STIS data reduction, see Tripp et al. 2005). In Fig. A18[ we show the COS and STIS normahzed 
profiles. While the STIS observations have a lower signal-to-noise level (typically S/N~ 10 per resolution element), the 
6.5 kms^^ FWHM resolution (similar to Keck HIRES Mg II observations) helps to characterize the velocity structure 
and saturation level. Unfortunately, the STIS cannot be used to determine the iVjji because the weak Lyman series 
transitions lie in a region of the spectrum where S/N< 3. Only C III, Si III, O VI, and H I are detected. The STIS 
Si III profile show thats a single component is present, which lines up very well with the weak H I transitions. The 
stronger H I transitions and C III and O VI show, however, additional components between about —150 and —25 
kms~^. For C III and O VI, these components are easily separated from the kms~^ component. For that absober, 
we were able to estimate the column densities in the narrow component at kms~^ for all the metal lines and weak 
H I transitions. The AOD estimates imply from COS and STIS Si III spectra yield 12.80 and 12.99 dex, implying that 
saturation effects are not negligible. Our adopted result for Si III is from a single Voigt component fit to the STIS 
profile, which yields log = 13.10 ± 0.09 and b = 9.9 ± 1.4 kms"""^. For C III, we can only derive a lower limit as the 
STIS profile is black. A single component COG fit to the weak H I lines (A916 to A923), the AOD measurement of 
same weak H I lines, and the continuum break method all give overlapping results for A'hi within la. We adopted a 
mean for A^hi of the results from these 3 methods. 

A Cloudy photoionization model with [X/H] < —1.65, \ogU < —2.6 reproduces simultaneously the column density 
of Si III and limits on Si II and Si IV. If Si III has a higher column density that derived from our profile fit, this would 
imply a higher U, but a lower metallicity (owing to the upper limit on Si II): therefore the adopted result implies 
a strict upper limit on the metallicity, but not on U. This model produces too much C III and N III. As C III is 
saturated, the discrepancy can be explained in terms of both saturation or nucleosynthesis effects. However, N III 
is no t detected, which implies that [N/Si] < —0.6, consistent with expected nucleosynthesis effect at low metallicity 
(e.g. Henry et al.||2000 ). Although the O VI directly associated with the LLS is weak and appears relatively narrow, 
the observed column density is several orders of magnitude predicted by the photoionization model; this absorber is 
therefore multiphase. 

From the results above and other studies at 2: < 1 (see Table IT]), we show in Fig. A19 the distribution of the central 
value of logU. This indicates a r elatively tight distributio n in U, with (log[/) — —3.3 ± 0.6 (median —3.2). There is 
an evolution of lo g U w ith z since Fumagalli et al. (2011a) repo rted that at 2; > 1.5, \ogU > —3. 

Finally, in Fig. A20 we provide a visual sumniary of Table Al for the A^hi, A^h, P/k, /, and T quantities as a 
function of the metallicity. There is no obvious trend between these quantities and the metallicity, except that the 
lowest values of A'h are observed in the high metallicity sample, which is expected (since the metallicity correction is 
smaller) . Importantly , the low temp eratures derived from the Cloudy simulations are confirmed by the observations 
(bottom panel of Fig. A20| and fjOl. 




Figure Al. Examples of Cloudy-predicted column densities as a function of ionization parameter U toward the same sightline PG1522+101 
for two different absorbers a,t z = 0.5185 {left panel) and z = 0.7289 {right panel). The gas is assumed to be photoionized by the background 
radiation field from quasars and galaxies, modeled as a uniform slab in thermal and ionization equilibrium. In both panels, the column 
density models are shown in the thinner lines, with the allowed values by the observations indicated by the thicker lines. The yellow region 
indicates the range of log U that fits best the observational constraints. Note that for the absorber on the right panel: i) a higher metallicity 
would predict too much Mg II for any U; ii) U is strictly a lower limit (higher U are allowed with lower metallicities); and ill) O VI is 
detected, but cannot be reproduced by photoionization. 



Table Al 

Cloudy Summary for LLS with 16.2 < logA^ni < 19 



Name 


z 


log A^H I 


logl/ 


[X/H] 


logA^H 


T 


lognn 


P/k 


/ 






[cm-2] 






[cm-2] 


(10" K) 


[cm"-^] 


(K cm~ 


3) (kpc) 


J0943+0531 


0.3544 


16.11 


< -3.2 


-1.3 


< 18.3 


< 1.5 


> -2.3 


> 79 


1.3 


J1419+4207 


0.4256 


16.17 


-2.9 


-1.4 


18.7 


1.7 


-2.5 


54 


0.5 


J1435+3604 


0.3878 


16.18 


> -3.0 


< -1.4 


> 18.6 


> 1.6 


< -2.4 


< 59 


0.3 


PG1116+215 


0.1385 


16.20 


-2.5 


-0.5 


19.1 


1.6 


-4.0 


2 


> 41 


PG1522+101 


0.5185 


16.22 


-3.6 


-0.4 


17.8 


1.1 


-1.9 


121 


0.02 


SBS1122+594 


0.5574 


16.24 


-3.1 


-1.0 


18.5 


1.5 


-2.2 


104 


0.2 


HE0439-5254 


0.6153 


16.28 


-2.7 


-0.3 


18.9 


1.2 


-2.5 


38 


0.8 


PG1338-I-416 


0.3488 


16.30 


-3.2 


-0.6 


18.5 


1.3 


-2.3 


68 


0.2 


J1419+4207 


0.5346 


16.34 


-3.9 


-0.2 


17.7 


0.7 


-1.4 


308 


0.004 


PG1407+265 


0.6828 


16.38 


-2.0 


-1.8 


20.0 


2.7 


-3.1 


19 


41 


PG1216-I-069 


0.2823 


16.40 


< -2.6 


< -1.7 


< 19.3 


< 1.9 


> -3.0 


> 21 


6.1 


J1419+4207 


0.2889 


16.40 


-3.1 


-0.6 


18.7 


1.3 


-2.4 


46 


0.4 


PG1338+416 


0.6865 


16.45 


-3.9 


0.1 


17.7 


0.6 


-1.2 


318 


0.003 


PKS0405-123 


0.1672 


16.45 


-3.1 


0.1 


18.9 


0.8 


-2.7 


40 


1.3 


J1619+3342 


0.2694 


16.48 


-2.9 


-1.6 


19.1 


1.7 


-2.7 


36 


2.0 


PKS0637-752 


0.4685 


16.48 


-3.8 


-0.4 


18.0 


0.9 


-1.5 


261 


0.006 


HE0153-4520 


0.2261 


16.61 


-3.8 


-0.8 


19.4 


1.9 


-2.9 


56 


6.4 


J1435+3604 


0.3730 


16.65 


-3.5 


-1.9 


18.6 


1.4 


-1.9 


156 


0.1 


PG1522+101 


0.7292 


16.66 


> -3.2 


< -2.0 


> 18.9 


> 1.5 


< -1.9 


< 190 


0.2 


PHL1377 


0.7390 


16.72 


-2.9 


-1.5 


19.3 


1.7 


-2.2 


105 


> 1.0 


PKS0552-640 


0.3451 


16.90 


> -3.8 


< -1.5 


> 18.5 


> 1.3 


< -1.7 


< 263 


0.05 


PG1630+377 


0.2740 


16.98 


-2.8 


-1.7 


19.7 


1.8 


-2.8 


30 


> 10 


PHL1811 


0.0810 


17.98 


-4.0 


-0.2 


19.0 


0.6 


-2.5 


19 


1.0 


PKS0312-77 


0.2026 


18.22 


-3.2 


-0.6 


19.8 


1.1 


-2.5 


35 


6.4 


TON153 


0.6610 


18.30 


-3.2 


-1.7 


19.4 


1.0 


-1.9 


141 


0.6 



Note. — z and A^hi are derived directly from the observations. U, [X/H], A'^jji T, rtj], P/k = n^T, and / = Nn/nn are derived from 
Cloudy photoionization simulations constrained by the column densities of metal ions and H I determined from the observations. For the 
sightlines that are not in the new sample, the results were retrieved from the original papers (see references in Table [ij. 
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Figure A2. Normalized absorption lines as a function of velocity centered on the absorber at 2 = 0.5185 toward PG1522+101. The UV 
transitions are from COS G130M and G160M and Mg II AA2796, 2803 are Keck HIRES spectra. The blue portion in each profile shows 
the approximate velocity range of the absorption (for Lya, the entire velocity interval is highlighted, while for the other species and other 
weak H I transitions, generally on the main H I component is highlighted). The reader should refer to Table [2] and the Appendix for the 
exact integration velocity intervals. The vertical dashed lines mark the zero velocity. 
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Figure A3. Same as Fig. |A2| but for the absorber at z 



100 -100 
= 0.7289 toward PG1522+101. 
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Figure A4. Same as Fig.p2^ but for the absorber &t z = 0.7390 toward PHL1377; C IV AA1548, 1550 are STIS E230M spectra. 
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Figure A5. Same as Fig.U2^ but for the absorber at at ^ = 0.3488 toward PG1338-I-416. 
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Figure A6. Same as Fig.lA2] but for the absorber at = 0.6865 toward PG1338+416. 
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Figure A7. Same as Fig.lA2] but for the absorber at 2 = 0.6828 toward PG1407+265. 
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Figure A8. Same as Fig.R2^ but for the absorber at 2 = 0.3451 toward PKS0552-640. 
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Figure A9. Same as Fig.p2^ but for the absorber at 2 = 0.4685 toward PKS0637-752. 
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Figure AID. Same as Fig.K2] but for the absorber at 2 = 0.6153 toward HE0439-5254. 
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Figure A12. Same as Fig.U2] but for the absorber at z = 0.2889 toward J1419+4207. 



34 



Lehner et al. 



m 
CD 



O 

CD 

N 
• I — I 

?^ 

o 



J14'iy+42'0^ "z='0|.4256' 




1.5 
1.0 
0.5 

0.0 
1.5 

1.0 

0.5 

0.0 
1.5 

1.0 

0.5 

0.0 
1.5 

1.0 

0.5 

Lg-i,!,i,-^-37-?. 

-200 



OVI A 1,031 




Silll A1206 , I 

• | - r i -mr i fiiii 1 1 1 1 1 [r r rrrrm- r i nii 1 1 1 1 1 r r r rrrr i - | -mr i i- 



MsII A2796 , 



1.0 lyii^^W^^M^'i'^^ 

0.5^ ! 



Figure A13. Same as Fig 
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21 but for the absorber at 2 = 0.4256 toward J1419+4207. 
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Figure A14. Same as Fig.U2] but for the absorber at z = 0.5346 toward J1419+4207. 
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Figure A15. Same as Fig.p2l but for the absorber at z = 0.3730 toward J1435+3604 
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Figure A16. Same as Fig.lXS] but for the absorber at z = 0.3878 toward J1435+3604. 
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Figure A17. Same as Fig.U2] but for the absorber at z = 0.2694 toward J1619-f3342. 
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Figure A18. Same as Fig. |A2[ but for the absorber at z = 0.2823 toward 
normalized STIS E140M profiTes. 
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Figure A19. Distribution of tlie central value of the ionization parameter, log U, derived for the LLS from the Cloudy photoionization 
models. 
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Figure A20. H I column densities determined from the observations and physical quantities determined from the Cloudy photoionization 
simulations as a function of the metallicity derived with the same models. In the bottom panel, the orange symbols show the result from 
the analysis of the temperature derived from the 6-value of the individual components of Mg II observed at high resolution with Keck 
HIRES. 



